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We present results of a numerical calculation of lattice QCD with two degenerate flavors of dynam- 
ical quarks, identified with up and down quarks, and with a strange quark treated in the quenched 
approximation. The lattice action and simulation parameters are chosen with a view to carrying 
out an extrapolation to the continuum limit as well as chiral extrapolations in dynamical up and 
down quark masses. Gauge conflgurations are generated with a renormalization-group improved 
gauge action and a mean field improved clover quark action at three values of /3 = 6/(?^, corre- 
sponding to lattice spacings of a ~ 0.22, 0.16 and 0.11 fm, and four sea quark masses corresponding 
to mps/mv ~ 0.8, 0.75, 0.7 and 0.6. The sizes of lattice are chosen to be 12^ x 24, 16^ x 32 
and 24^ x 48 so that the physical spatial size is kept constant at La « 2.5 fm. Hadron masses, 
light quark masses and meson decay constants are measured at five valence quark masses corre- 
sponding to mps/mv ~ 0.8, 0.75, 0.7, 0.6 and 0.5. We also carry out complementary quenched 
simulations with the same improved actions. The quenched spectrum from this analysis agrees well 
in the continuum limit with the one of our earlier work using the standard action, quantitatively 
confirming the systematic deviation of the quenched spectrum from experiment. We find the two- 
flavor full QCD meson masses in the continuum limit to be much closer to experimental meson 
masses than those from quenched QCD. When using the K meson mass to fix the strange quark 
mass, the difference between quenched QCD and experiment of 2.61q g% for the K* meson mass 
and of 4.lt?:6% for ^he </> meson mass is reduced to 0.7^^ 7% and 1.3l2;5% in QCD, where the 
errors include estimates of systematic errors of the continuum extrapolation as well as statistical 
errors. Analyses of the J parameter yield a similar trend in that the quenched estimate in the 
continuum limit J = 0.375lg QQg increases to J = 0.440lg Q3} in two-fiavor full QCD, approaching 
the experimental value J ~ 0.48. We take these results as manifestations of sea quark effects in 
two-fiavor full QCD. For baryon masses full QCD values for strange baryons such as H and are 
in agreement with experiment, while they differ increasingly with decreasing strange quark content, 
resulting in a nucleon mass higher than experiment by 10% and a A mass by 13%. The pattern 
suggests finite size effects as a possible origin for this deviation. For light quark masses in the 
continuum limit we obtain m^(2 GeV) = iAAt^lt MeV and m^(2 GeV) = 88l g MeV (if-input) 

and m^^(2 GeV) — 90l^j MeV (^-input), which are reduced by about 25% compared to the values 
in quenched QCD. We also present results for decay constants where large scaling violations ob- 
struct a continuum extrapolation. Need for a non-perturbative estimate of renormalization factors 
is discussed. 



* address tiU 31 August, 2000 

^present address : Department of Physics, University of Wales, Swansea SA2 8PP, U.K. 
■'■present address : Physics Department, Columbia University, New York, NY 10027, U. S. A. 
^present address : CERN, Theory Division, CH-1211 Geneva 23, Switzerland 

**present address : Department of Biochemistry and Molecular Biology, University College London, London, England, U.K. 



1 



I. INTRODUCTION 



The mass spectrum of hadrons represents a fundamental manifestation of the long-distance dynamics of quarks and 
gluons governed by QCD. Non-perturbative calculations through numerical simulations on a space-time lattice 0] 
provide a method to obtain this quantity from the QCD Lagrangian without approximations. Such calculations also 
lead to a determination of the light quark masses , which are fundamental constants of Nature and yet not directly 
measurable in experiments. These reasons underlie the large number of attempts toward the hadron spectrum carried 
out since the pioneering studies of Ref. |^ . 

Most of these calculations employed the quenched approximation of ignoring sea quarks, since dynamical quark 
simulations incorporating their effects place quite severe demands on computational resources. Significant advance 
has been made over the years within this approximation. In particular, Weingarten and collaborators [Q made a 
pioneering attempt toward a precision calculation of the spectrum in the continuum limit through control of all 
systematic errors other than quenching within a single set of simulations. 

This approach was pushed further in Ref. Q where the precision of the calculation reached the level of a few percent 
for hadron masses. Scrutinized with this accuracy, the quenched hadron spectrum shows a clear and systematic 
deviation from experiment; when one uses tt, p and K meson masses as input to fix the physical scale and light quark 
masses, the K* — K hyperfine splitting is too small by about 10% compared to the experimental value, the octet 
baryon masses are systematically lower, and the decuplet baryon mass splitting is smaller than experiment by about 
30%. 

Clearly further progress in lattice calculations of the hadron mass spectrum requires a departure from the quenched 
approximation. In fact simulations of full QCD with dynamical quarks have a long history Joj-p^, leading up to the 
recent efforts of Refs. p^-[T^]. In contrast to quenched simulations, however, no attempt to control all of the systematic 
errors within a single set of simulations has been made so far. Except for the work of the MILC Collaboration [ p^ , 
employing the Kogut-Susskind quark action, previous calculations have been restricted to a few quark masses within 
a small range and/or a single value of the lattice spacing. Furthermore, until recently, statistics have been rather 
limited due to the limitation of available computing power. 

In the present work, we wish to advance an attempt toward simulations of full QCD which includes extrapolations to 
the chiral limit of light quark masses and the continuum limit of zero lattice spacing. This is an endeavor demanding 
considerable computing resources, which we hope to meet with the use of the CP-PACS parallel computer with a peak 
speed of 614 GFLOPS developed at the University of Tsukuba [ p9|j20| ]. We explore, as a first step toward a realistic 
simulation of QCD, the case of dynamical up and down quarks, which are assumed degenerate, treating the strange 
quark in the quenched approximation. Preliminary results of the present work have been reported previously pl|] . 

A crucial computational issue in this attempt is how one copes with the large amount of computation necessary 
in full QCD, and still covers a range of lattice spacings required for the continuum extrapolation. We deal with this 
problem with the use of improved lattice actions, which are designed to reduce scaling violations, and hence should 
allow a continuum extrapolation from coarse lattice spacings. 

In Ref. we have carried out a preparatory study on the efficiency of improved actions in full QCD. Based on 
the results from this study we employ a renormalization group improved action [ p3| for the gauge field and a mean 
field improved Sheikholeslami-Wohlert clover action for the quark field. With these actions, hadron masses show 
reasonable scaling behavior and the static quark potential good rotational symmetry, at a coarse lattice spacing of 
a « 0.2 fm, as compared to the range a<0.1 fm needed for the standard plaquette and Wilson quark actions. This 
leads us to make a continuum extrapolation from the range of lattice spacings a « 0.2-0.1 fm. 

Previous studies of finite size effects (see, e.g., Refs. P JlT| , p^ ) indicate that physical lattice sizes larger than 
La « 2.5-3.0 fm are required to avoid size-dependent errors in hadron masses. Compromising on a lattice of physical 
size La ~ 2.5 fm leads to a 12"^ x 24 lattice at a ~ 0.2 fm, and 24^^ x 48 at a w 0.1 fm. Estimates of CPU time 
obtained in our preparatory study [p2[ show that simulations on such a set of lattices are feasible with the full use of 
the CP-PACS computer. 

Since we employ the quenched approximation for strange quark, the calculation of the strange spectrum requires 
the introduction of a valence strange quark which only appears in hadron propagators. We generalize this treatment 
and analyze hadron masses as functions of valence and sea quark masses regarded as independent variables. The 
benefit of this approach is that it gives us better control over the whole spectrum (strange and non-strange) and its 
cross-over from quenched to full QCD when the mass of the underlying sea quark is decreased. 

There are a number of physics issues we wish to explore in our full QCD simulation. An important question is 
whether effects of dynamical quarks can be seen in the light hadron spectrum. In particular we wish to examine if and 
to what extent the deviation of the quenched spectrum from experiment established in our extensive study with the 
standard plaquette and Wilson quark actions ||] can be explained as effects of sea quarks. Answering this question 
requires a detailed comparison with hadron masses in quenched QCD for which we use results of Ref. We also 
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carry out a set of new quenched simulations with the same RG-improved gluon action and the clover quark action as 
employed in the simulation of full QCD in order to make a point-to-point comparison of full and quenched QCD at 
the same range of lattice spacings. 

Another question concerns light quark masses. Quenched calculations of light quark masses have made considerable 
progress in recent years p8|J^] . It has become clear that the quenched estimate for the strange quark mass 
extrapolated to the continuum limit suffers from a large systematic uncertainty of order 20% depending on the choice 
of hadron mass for input, e.g., K meson mass or (p meson mass. This is a reflection of the systematic deviation of the 
quenched spectrum from experiment. It is an important issue to examine how dynamical quarks affect light quark 
masses and resolve the systematic uncertainty of strange quark mass. A recent attempt at a full QCD determination 
of light quark masses [g9[ was restricted to a single lattice spacing. We extracted light quark masses through analyses 
of hadron mass data obtained in the spectrum calculation. The main findings of our light quark mass calculation 
have been presented in Ref. |3^. We give here a more detailed report of the analysis and results. 

Full QCD configurations generated in this work can be used to calculate a large variety of physical quantities 
and examined for sea quark effects. We have already pursued calculations of several quantities. Among these, 
the flavor singlet meson spectrum and its relation with topology and U{\) anomaly is of particular interest from 
the theoretical viewpoint, and preliminary results have been published in Ref. Other calculations concern the 

prediction of hadronic matrix elements important for phenomenological analyses of the Standard Model. Results have 
been published for heavy quark quantities such as B and D meson decay constants |^,^^ as well as bottomonium 
spectra [Q. A report of the analysis of the light pseudoscalar and vector meson decay constants is included in this 
article. 

The outline of this paper is as follows. We first describe details of the lattice action, the choice of simulation 
parameters and the algorithm for configuration generation in Sec. O. Measurements of hadron masses, the static 



quark potential and a discussion of autocorrelations are presented in Sec. [II. In Sec. IV we discuss the procedure of 
chiral extrapolation. Section ^ contains the main results for the full QCD light hadron spectrum. In Sec. ^ we then 
turn to a presentation of quenched QC D sim ulations with improved actions. This sets the stage for a discu ssion of sea 



quark effects which is contained in Sec. VII. Calculations of light quark masses are presented in Sec. VII]| . Section 



contains a discussion of decay constants. Finally, we present our conclusions in Sec. ^ 



IX 



II. SIMULATION 

A. Choice of improved lattice action 

Based on our preparatory study in Ref. we choose improved gauge and quark actions for full QCD configuration 
generation. The improved gluon action has the form 



(1) 

The coefficient ci = —0.331 of the 1x2 Wilson loop W^^"^ is fixed by an approximate renormalization group 
analysis and cq = 1 — 8ci = 3.648 of the 1x1 Wilson loop by the normalization condition, which defines the 
bare coupling (3 = 6/5^. From the point of view of Symanzik improvement the leading scaling violation of this action 
is O(a^), the same as for the standard action. 

For the quark part we use the clover quark action p4| ] defined by 

Sq^^qxDx.yqy, (2) 

'x,y 

Dx,y ^ Sxy - {{1 - lfj)Ux,fj,6x+ij.,y + (1 + lli)Ul^^5x'y+ii] ~ SxyCSW^'^ Crfj,iyFf,iy, (3) 

where k is the usual hopping parameter and F^j^u the standard lattice discretization of the field strength. 
For the clover coefficient csw we adopt a mean field improved choice defined by 

csw = (I^ixi)-'/' = (1 - Q.MUrT'" ' (4) 

where for the plaquette W^^^ the value calculated in one- loop perturbation theory |^ is substituted. This choice 
is based on our observation in Ref. |22| that the one-loop calculation reproduces the measured values well. Indeed, 
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an inspection of Table XXIV in Appendix |^ shows that ^ ^ in the simulations agree with one- loop values with a 
difference of at most 8%. The agreement for csw is not fortuitous; the one-loop value for the RG gauge action (0), 
which was calculated after the present work was started, equals csw = 1 + 0.678//3, which differs from our choice 
Csw = 1 + 0.631//3 + . . . only by a few per cent. We do not employ the measured plaquette for the clover coefficient, 
as prescribed by the usual mean field approximation, which would have required a time-consuming self-consistent 
tuning. The leading scaling violation with our choice of csw is 0{g^a). 



B. Simulation parameters 

The target of this work is a calculation of the two-flavor QCD light hadron spectrum in the continuum limit and at 
physical quark masses. For this purpose we carry out simulations at three lattice spacings in the range a « 0.2-0.1 fm 
for continuum extrapolation, and at four sea quark masses corresponding to mT^/mp « 0.8-0.6 for chiral extrapolation 
for each lattice spacing. The simulation parameters are given in Table |. 

We employ three lattices of size 12^ x 24, 16"^ x 32 and 24^^ x 48 for our runs. The coupling constants (3 = 1.8, 1.95 and 
2.1 are chosen so that the physical lattice size remains approximately constant at La « 2.5 fm. The resulting lattice 
spacings determined from the p meson mass equal a = 0.2150(22), 0.1555(17) and 0.1076(13) fm or a^^ = 0.9177(92), 
1.269(14) and 1.834(22) GeV. 

We have also performed an initial run at /3 = 2.2 on a 24^ x 48 lattice for which the lattice spacing turned out 
to be a = 0.087 fm. The physical lattice size La = 2.08 fm is significantly smaller than the other three lattices. In 
order to avoid different magnitude of possible finite size effects, we do not include data from this run when we make 
extrapolations to the continuum limit. They will be included in figures and tables for completeness, however. 

We carry out hadron mass analyses distinguishing the sea and valence quark hopping parameters Ksoa and K^ai ■ At 
each value of /3, configurations are generated at four sea quark hopping parameters Kgoa such that the mass ratio of 
pseudoscalar to vector mesons made of sea quarks takes mps/my ~ 0.8, 0.75, 0.7 and 0.6. At each sea quark mass, 
hadron propagators are measured for five valence hopping parameters Kvai with approximate ratios of mps/niy « 0.8, 
0.75, 0.7, 0.6 and 0.5. The four heavier Kvai coincide with those chosen for sea quarks. 

A schematic representation of our choice on the (l/^sca, l/^vai) plane is shown in Fig. ^ The physical point is 
characterized by l/Kgoa = l/'^vai = ^/i^ud for degenerate up and down quarks, and l/Kgoa = ^/I'^ud and l/^vai — 
l/^strango for straugc quark, i.e., lying in the shaded region on the left hand side of the diagram. The additional 
points with l/nvai = V5 in the bottom part of the diagram are not directly needed in exploring the physical region. 



As we will see in Sec. IV, they help in the description of hadron masses as a combined function of sea and valence 
quark masses and are therefore indirectly useful for the extrapolation to physical points. Including them also keeps 
the possibility open for a future extension of the present work towards the chiral limit by adding the fifth sea quark 
and completing the grid of Fig. |l]. 

Our choice of hopping parameters enables us to obtain the full strange and non-strange hadron spectrum in a sea 
of degenerate up and down quarks. If we denote with S a valence quark with Kvai = Ksea and with V a valence quark 
with Kvai ^ Ksca, we obtain mesons of the form SS, SV and VV and baryons of the form SSS, SSV, SVV and VVV. 



C. Configuration generation 

Configurations are generated for two fiavors of degenerate quarks with the Hybrid Monte Carlo (HMC) algorithm. 
In Table || we give details of the parameters and statistics of the runs. At the main coupling constants (3 — 1.8-2.1, 
runs are made with a length of 4000-7000 HMC unit-trajectories per sea quark mass. The add itional runs at /3 = 2.2 



are stopped at 2000 HMC trajectories per sea quark mass for the reason described in Sec.|IB. 

To speed up the calculation we have implemented several improvements in our code. For the inversion of the quark 
matrix during the HMC update we use the even/odd preconditioned BiCGStab algorithm [|6|. Test runs confirmed 
that the performance of this algorithm is better than that of the MR algorithm and that the advantage increases 
toward lighter quark masses ]37| ]. In a test run at m,r/wp ~ 0.7 we observed a 43% gain in computer time for the 
same accuracy of inversion compared to the MR algorithm. 

As a stopping condition for the inversion of the equation D{k)G ~ B during the fermionic force evaluation we use 
the criterion 

\\DG - BW^ < stop, (5) 

with values of stop given in Table ^ where we also give the number of iterations necessary for the inversion. For the 
evaluation of the Hamiltonian we use a stricter stopping condition which is smaller by a factor of 10^ than the one 
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used for the force evaluation. With these stopping conditions, the Hamiltonian is evaluated with a relative error of 
less than 10^^°. We have also checked that the reversibility over trajectories of unit length is satisfied to a relative 
level better than 10^^ for the gluon link variables. 

Another improvement concerns the scheme for the integration of molecular dynamics equations. For our runs we 
have used the following three schemes. 

a) the standard leap-frog integration scheme: The operator to evolve gauge fields and conjugate momenta by a step 
At in fictitious time can be written in the form 



Tp 



(6) 



where the operator Tp{At) — exp(AT Pidi) moves the gauge field U hy a, step At, whereas the operator Tq(At) = 
exp(— At diS{U, ^)d/dpi) moves the conjugate momenta p by a step At. The leap-frog integrator has an error of 
O(At^) for a single step and of O(At^) for a unit-trajectory. 

b) an improved scheme: The discretization error of the leap-frog integration scheme can be reduced by using an 
improved scheme. The simplest improvement has the form. 



Tp -At Tc 



^)tp((1-5)At)T,(^ 



Tp t^At 



(7) 



This scheme has errors of the same order as the standard leap-frog scheme but the main contribution to the error is 
removed by the choice 6 = (3 — •\/3)/3 Test runs have shown that At can be taken a factor 3 larger than for 

leap-frog without losing the acceptance rate for the heaviest sea quark. This leads to a gain of about 30% in computer 
time. The gain, however, decreases toward lighter quark masses, and the computer time required for the improved 
scheme at the lightest quark mass is roughly the same as for the standard leap-frog scheme (see parameters of the 
run at /3 = 1.8 and k = 0.1464 in Table |l|). 

c) Sexton- Weingarten scheme Q: In this scheme the evolution with the gauge field force diSg{U) is made with 
an n times smaller time step than that with the fermionic force diSf{U, $) according to 



At 
'2n 



1^2 (At 



Ti 



At 
"2n 



(8) 



where 



Ti (At) = Tp ( -At 



exp —A? 



^d,Sg{U)d/dp,\Tp 



-An 



T2 (At) = exp (^-ATY,^^Sf{U,<P)^/^p, 



(9) 



(10) 



We have implemented a scheme for which both Eq. (^) and Eq. (^) are improved as in Eq. (^ . For n = 2 the time step 
At can be chosen 10% larger than in scheme b) while maintaining a similar acceptance. However, this improvement 
is offset by an increase of a factor 4 in the number of operations for the gauge field force. This leads to an increase of 
30% in the total number of operations at ^ = 1.8, k — 0.1445. Hence the performance of scheme c) is similar to the 
leap-frog scheme, as can be seen in Table ||. 

The scheme employed for each run is listed in Table |l]. After some trials on the smaller lattices (12^^ and 16'^) we 
found the scheme b) to be most practical, and we used it for all the runs on the larger 24^ lattices. The step size At 
for molecular dynamics has been chosen so that the acceptance ratio turns out to be 70-80%. 

Light hadron propagators are measured simultaneously with the configuration generation with a separation of 5 
HMC trajectories. The number of measurements is given in Table ||. We stored configurations with a separation of 
10 HMC trajectories (at (3 — 1.8 and 1.95) or 5 HMC trajectories (at (3 — 2.1 and 2.2) on tapes for later measurement 
of other observables such as the topological charge and flavor singlet meson mass 31 1 , quarkonium spectra and 
the B meson decay constant [ ^|33| . 

In the last column of Table y , we list the number of configurations removed by hand because of the occurrence of 
exceptional propagators. We did not encounter exceptional configurations in full QCD where Kvai = ^soa- However, 
strange behavior of propagators did occur for the lightest valence quark mass for some configurations. We have 
removed all the propagators obtained on such configurations in order to allow a jack-knife error analysis. 

Our criterion for removal of a configuration is a deviation of hadron propagator by more than 10 standard deviations 
from the ensemble average for at least one channel and at least one timeslice. The fraction of removed configurations 
drops from 1.2% at (3 — 1.8 to 0.1% at /3 = 2.1. No configurations needed to be removed a.t (3 = 2.2. 
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D. Coding and runs on the CP-PACS computer 



We have spent much effort in optimizing the double precision codes for configuration generation on the CP-PACS 
computer as described in Ref. |3^]. Actual runs took advantage of the partitioning capability of the CP-PACS, using 
64 PU (processing units), 256 PU and 512 PU for the lattice size 12^ x 24, 16^ x 32 and 24^ x 48, and executing runs at 
different values of Ksea at the same time. For some of the runs at smaller quark masses, which need longer execution 
times, we made two or more independent parallel runs which are combined for the purposes of measurements. 

The CPU time needed per trajectory is listed in Table ||. Converted to the number of days with the full use of the 
CP-PACS computer, the configuration generation took 10 days for (3 — 1.8 on a 12^^ x 24 lattice, 40 days at /3 = 1.95 
on a 16"^ X 32 lattice, 186 days at /3 = 2.1 on a 24^^ x 48 lattice and 82 days on the same size lattice at /3 = 2.2. Adding 
3+12+46+23 days for measurements of observables and 1+3+6+3 days for I/O loss, the entire CPU time spent for 
the simulations equals 415 days of the full operation of the CP-PACS computer. 



III. MEASUREMENTS 



A. Hadron masses 



We employ meson operators defined by 

M{<^{n)^7„rA9n, (11) 

where / and g are quark fields with flavor indices / and g, and Ta represents one of the 16 spin matrices Ta = I, 75, 
*7/i75: 7p and «[7^,7iy]/2 of the Dirac algebra. Using these operators, meson propagators are calculated as 

{MA{n)MAm- (12) 

For the operator of octet baryons with spin J = 1/2 we use the definition 

Oi^\n) = e^'^if^^Cj.gtW^^, (13) 

where a, 5, c are color indices, C = 7472 is the charge conjugation matrix and a = 1, 2 represents the z-component of 
the spin = ±1/2. To distinguish E and A-like octet baryons we antisymmetrize fiavor indices, written symbolically 
as 

^^jmi+Ml, (14) 

v2 

[fh]g-[gh]f-2[fg]h ^ ^^^^ 

where [fg] fg - gf. 

The operator of decuplet baryons with spin J = 3/2 is given by 

Di'aH - e'^'^if^'^Cj^g'jhl^. (16) 

Writing out the spin structure (^, a) explicitly, we employ operators for the four z-components of the spin Jz — 
±3/2, ±1/2 defined as 

Dy2 = e'^^^{f'^CV+g')hl, (17) 
Di,2 = e^'-'lif^CT^g'Wi - {f^CV+g'')hl]/:S, (18) 
^-1/2 = e^'^'^lif^CT^g'^M - {f''CT^g'^)hl]/i, (19) 
D-i/2 = e"''^(/^"Cr_g^)/i^, (20) 

where r± = (71 + i72)/2 and Fq = 73. 

Using operators defined as above, we calculate 8 baryon propagators given by 

(I]„(n)E„(0)), a =1,2, (21) 
(A„(n)A,(0)), a = 1,2, (22) 
{Ds{n)Dsm. 5 = 3/2,1/2,-1/2,-3/2, (23) 
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together with 8 antibaryon propagators similarly defined. 

We average zero momentum hadron propagators over three polarization states for vector mesons, two spin states 
for octet baryons and four spin states for decuplet baryons. We also average the propagators for the particles with 
the ones for the corresponding antiparticles. 

For each configuration quark propagators are calculated with a point source and a smeared source. For the smeared 
source we fix the gauge configuration to the Coulomb gauge and use an exponential smearing function ■(/'('') = 
Aexp{—Br) for r > with V^fO) — 1. We chose A and B based on experiences from previous quenched measurements 
of the pion wave function [po[ and from our preparatory full QCD study |p2| and readjusted them by hand so that 
hadron effective masses reach a plateau as soon as possible on average. The values of A and B are given in Table III. 

In Figs. I^-IJ we show examples of effective mass plots for hadron propagators with degenerate valence quarks equal 
to the sea quark. Effective masses from hadron propagators where all the quark propagators have been calculated 
with smeared sources have the smallest statistical errors and exhibit good plateaux starting at smaller values of t than 
those containing point sources. We therefore use smeared propagators for hadron mass fits. 

Fit ranges [tmin,imax] are determined by inspecting effective mass plots. As a general guideline, we choose <min 
to be in the same range for given particle and gauge coupling. The approach to a plateau changes with the quark 
mass and this allows for a variation of tmin- To be confident that contributions of excited states die out at we 
also consult effective masses from propagators with point and mixed sources. The upper end of the fit range, imax, is 
chosen to extend as far as the effective mass exhibits a plateau and the signal is not lost in the noise. 

Hadron masses are derived from correlated fits to propagators with correlations among different time slices taken 
into account. We assume a single hyperbolic cosine for mesons and a single exponentia l for b aryons. With a statistics of 
4000-7000 HMC trajectories (corresponding to 80-140 binned configurations, see Sec. HID) for hadron propagators at 
/3 = 1.8, 1.95 and 2.1, the covariance matrix is determined well. Typically, the errors of eigenvalues of the covariance 
matrix are around 15%, and fits have a x^/-^df around 1 and at most 3. For /? = 2.2, however, where fewer 
configurations are available, eigenvalues of the covariance matrix have typical errors of 30%, and the correlated fits 
are less stable. For all the cases we also made uncorrelated fits and checked that masses are consistent within error 
bars. 

Errors in hadron masses and in x^/^Ydf are estimated with the jack-knife procedure with a bin size of 10 configu- 



rations or 50 HMC trajectories. A discussion of the choice of this bin size will follow in Sec. |III D . 

Resulting hadron masses are collected in Appendix There and in the following, lower case symbols are used for 
observables in lattice units, for which the lattice spacing a is not explicitly written. 



B. Quark mass 

Another quantity which can be obtained from meson correlation functions is the quark mass based on the axial 
vector Ward identity ( AWI) 1^,^ . It is defined from matrix elements of the pseudoscalar density P and the fourth 
component of the axial vector current A^ by the expression 

_ (oiv4AriPS) 

" 2(0|F|PS) ' ^ ' 

where we employ the improved axial vector current A™^ = A4 + CAd^P with ca calculated in one-loop perturbation 
theory and representing the symmetric lattice derivative (see Appendix |c|) . 

In practice we extract the AWI quark mass from single-exponential fits to meson correlators. For the analysis of 
pseudoscalar masses we assume the form 

(P(i)P(O)) = Cp [cxp(-mpsi) + exp(-mps(it - t))] , (25) 

which has already been described above. Keeping the value of mpg obtained from this fit, we make an additional fit 
to the correlator 

(4'"P(i)P(0)) - Ca [exp(-mpst) - exp(-mps(Lt ~ t)] , (26) 
where Ca is the only fit parameter. The AWI bare quark mass before renormalization is then obtained through 

AWI _ ™psCa 



2Cp 

Results for are given in Appendix 



(27) 
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C. static quark potential 



We measure the temporal Wilson loops applying the smearing procedure of Ref. |Q| . The number of smearing steps 
is fixed to 2, 4 and 10 on 12'^ x 24, 16^ x 32, and 24^^ x 48 lattices, respectively, which we find sufficient to ensure a 
good overlap of Wilson loops onto the ground state. The static quark potential V{r) is determined from a correlated 
fit of the form 

W{r,t) ^C{r)exp{-V{r)t) . (28) 
As shown in Fig. || noise dominates the signal when the temporal size of W{r, t) exceeds t w 0.9 fm. We therefore 



take fit ranges, listed in Table IV, which approximately correspond to t ^ 0.45-0.90 fm at /3 = 1.8, 1.95 and 2.1. At 
f3 = 2.2, we use the same fit ranges as those taken at /3 = 2.1. 

A typical result for V{r) is plotted in Fig. ^. Since we do not observe signs of string breaking, we parametrize V{r) 
in the form, 

V{r) = Vb - - + err. (29) 
r 

The lattice correction to the Coulomb term calculated from one lattice gluon exchange diagram nm is not included 
since breaking of rotational symmetry is found to be small with the improved actions we employ [ p2[ . 
The Sommer scale tq is defined through |H] 



dV{r) 



Using the fit parameters in Eq. (ES), tq can be obtained from 



1.65. (30) 



; 1.65 — a , 
^o = \ ■ (31) 



We fit potential data to Eq. (BS) and determine ro for several fitting ranges lying in the interval [Rmin,Rj 



maxj 



Values of i?min and i?max are listed in Table |V[ We take the average of fit results as central values for Vb, a, cr and 
To, and use the standard deviation as an estimate of the systematic error. Results of a and ro are summarized in 
Table 0. 

D. Autocorrelations 

The autocorrelation function of a time series of a variable / is defined as 

where the unnormalized autocorrelation function is given by 

Cfit) = - (33) 

in 



The quantity relevant for the determination of the statistical error of / is the integrated autocorrelation time rl,"*. 



defined as 



-. CXD -. OO 

^)"*-o E p/W-o+E^/W- (34) 



2 ^ ' 2 

t— — oo t—1 



The naive error estimate is smaller than the true error by a factor of J 2t|"* . In numerical estimations of rj"* , the 

sum in Eq. (jsj) has to be cut off. It has been found to be practical to calculate the sum self-consistently up to 
t « (4-10)rj:"*. A convenient quantity for this purpose is the cumulative autocorrelation time 



1 * 

W = ^+Ep/(^)' (35) 



2 .=1 



which should run into a plateau for r^"™(t) w i/4-i/lO. 

We calculate autocorrelation times for three different quantities: 
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1. The gauge action cqW + ciW . Measurements are made after every HMC trajectory. 

2. The number of iterations iVinv for the inversion of the Dirac matrix during the HMC update. Since this quantity 
is governed by the ratio of the largest to the smaUest eigenvalue of the Dirac matrix, it is expected to be the 
quantity which takes the longest simulation time to decorrelate. Measurements are made during every HMC 
trajectory. 

3. The effective pion mass mTr.cff measured at the onset of a plateau. Measurements are made only after every 5**^ 
HMC trajectory. 

Two examples for autocorrelation function and cumulative autocorrelation time are shown in Fig. |^. The cumulative 
autocorrelation time shows a plateau around the expected region from which we estimate the integrated autocorrelation 



times Tf given in Table VI 



Values of r"^* are generally below 10 HMC trajectories for the runs at /? < 2.10. These numbers are significantly 
lower than initial estimates for the HMC algorithm ^ and also lower than estimates from recent simulations with 
the Wilson or clover fermion action 0,0. A possible reason might be coarser lattice spacings of our simulations 
compared to the studies mentioned above. It has also been noticed in Ref. ||l^ that autocorrelations appear to be 
weaker on larger lattices. Our lattice sizes in physical units are considerably larger than the ones in Refs. |]l7| , ^ . 

Another point of interest is the size of increase of the autocorrelation time with decreasing sea quark mass. For the 
gauge action and for Ninv the autocorrelation time grows by about a factor of two in the range of simulated sea quark 
masses, whereas for the effective pion mass the situation is less clear. These observations are roughly consistent with 
the findings in Refs. 

A practical way to take into account autocorrelations in error analyses is to use the binning method. In Fig. ^ 
we show the increase of the relative error of the pion mass as a function of the bin size. The plotted error bars are 
determined by a jack-knife on jack-knife method. For this plot we have used uncorrelated fits to the pion propagator, 
since for larger bin sizes the number of configurations would not be large enough to reliably determine the covariance 

matrix for correlated fits. We observe that the error rises to a plateau which is about a factor ^2t™* larger than the 

naive error obtained with a unit bin size. From these and similar figures at other simulation parameters we find that 
a bin size of 10 configurations, equivalent to 50 HMC trajectories, covers all the autocorrelations we have examined 
while leaving sufficient number of bins to allow correlated fits. We therefore employ this bin size in all error analyses. 



IV. CHIRAL EXTRAPOLATIONS 



The calculation of the physical hadron spectrum requires an extrapolation from simulated quark masses to the 
physical point. In order to make these extrapolations we have to fit hadron masses to a functional form chosen to 

(i) 

express their chiral behavior. Hadron masses are functions of Kgea and k^^i, where i = 1,2,... labels valence quarks. 
We take this into account by performing combined fits to all measured masses of a given channel. 

The hopping parameter is not the only choice for the basic variable in these fits. Pseudoscalar meson masses can 
be employed as well for vector mesons and baryons. This has the advantage that only measured hadron masses are 
involved, and we employ this way of parametrizing vector meson and baryon masses. Pseudoscalar meson masses 
themselves, however, have to be expressed in terms of quark masses in order to fix the physical point in terms of 
quark masses. 



A. Pseudoscalar mesons 



Let us recall that the definition of quark mass suggested by a Ward identity for vector currents (VWI) has the form 

where Kc is the critical hopping parameter at which pseudoscalar meson mass vanishes. For a combined fit of 
pseudoscalar meson masses in terms of this "VWI" quark mass, we define sea and valence quark masses through 

^ywi ^ 1 fj_ _ 1\ ^ (37) 
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where k'^^^ denote for i = 1, 2 the hopping parameters of the valence quark and antiquark which make the meson. In 
the leading order of chiral perturbation theory the masses squared of pseudoscalar mesons are linear functions of the 
average quark mass. We therefore define an average valence quark mass through 

Figure ^ shows pseudoscalar meson masses as functions of 1/Kvai- We observe that partially quenched data {i.e., 
VV and SV) lie along clearly distinct lines when the hopping parameter of sea quark Ksc^ is varied. Each of the 
partially quenched data are close to linear, but their slope shows a variation with Ksca,- As illustrated in the inlaid 
figures we also see that the VV and SV masses lie along slightly different lines, which means that masses depend on 
the individual valence quark masses rn^^^j^ in addition to their average. 

These features of pseudoscalar meson mass data lead us to adopt a fit ansatz which consists of general linear and 
quadratic terms in the valence quark mass and in the sea quark mass given by, 



tpg 1 f„soa, '^val' '^val 



'*val''^valj — "s^^soa +"^771^^^ + Cg [m^^^ ) + Cy [m^^i ) 



_VWI VWI , VWI ™VWI /.(.x 

^"^soa "^val + COT"^val(l)'^val(2)■ ) 



Figure || shows the fit with solid lines for the SS channel and with dashed (SV) or dot-dashed (VV) lines for partially 
quenched data. The lines follow the data well. We employ uncorrelated fits for chiral extrapolations even though 
data with common Ksea are expected to be correlated. Obtained values of x^/A^df can therefore only be considered 
as rough guidelines to judge the quality of fits. Except for (3 — 1.8 where x^/iVoF = 4, we obtain values which are 



smaller than 1. Fit parameters Kc, b's for linear terms and c's for quadratic terms and x^/-^df are given in Table VI] 



A different definition of quark mass suggested by a Ward identity for axial vector currents is given by Eq. (E3) . Since 

(i)T 
van 



this is a measured quantity derived from meson propagators it depends on three hopping parameters, K^^J^^ii — 1,2 



of the valence quark and antiquark, and Ksca of the sea quark. We define 

(«sca;4142i), (41) 

(42) 

^sea ~ ™ (/^sca; '^scaj '^sca) ■ (43) 

Pseudoscalar meson masses are expressed in terms of these quantities with the quadratic ansatz, 

,..(1) J^)\ -V^AWI , / Cm^Wn2 / ^AWI AWI (..^ 



m^WI _ ^AWI 



™AWI 
"T'val 2 



'PS 



In contrast to Eq. (^) monomial terms in the sea quark mass are absent since pseudoscalar masses vanish in the 
chiral limit ra^^ = for each value of the sea quark mass. Data of different degeneracies lie on common l ines and 



therefore we have dropped the term with individual "^^^Jy Fit parameters and x^/-^df are given in Table VIl 



B. Vector mesons 



Vector meson masses are fit in terms of measured pseudoscalar meson masses. We define 

/i, = m|g (^Ksca; 4al' '^val) ' (45) 
Alval = ^ (/il +Ai2) , (46) 
A'sca — ™PS ('*scai ^scaj '^soa) ■ (47) 

Vector meson masses as functions of /ivai are shown in Fig. The general feature of data is similar to the one for 
pseudoscalar mesons. We find, however, that the lines for VV and SV are indistinguishable. Hence, vector meson 
masses do not require terms in individual /i^'s. We therefore take a quadratic function in /^soa and /^vai of the form 

mv («.ca; 4al: 4al) = + + S^A^val + C^A^^a + C^Mval + /^scaMval • (48) 
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Fit li nes de scribe data well as shown in Fig. nol and x^/^df is at most 1.4. Fit parameters and x'^/^'df are given in 



Table VIII 



Chiral perturbation theory predicts | |48| that the first correction to the linear term in /i has a non-analytic 3/2 
power of /i. In order to examine if data show evidence for such a dependence, we attempt a fit of the form 

/ (1) (2)\ ,v I T-,V I rjV I j-^V 3/2 , j-^V 3/2 „v 1/2 

my (^Ksoa; <al' <al j ^ + Msoa + ^val + Msca + ^.J^y + ^'^^MsoaMval ■ (49) 

1 /2 . . . 1/2 

The cross-term of the form /isoaMvai gi'^^s rise to a term proportional to /i^l^j for the partially quenched case where 
/isoa is kept constant. This is similar to quenched QCD. Terms proportional to /ig/a are expected to be absent |49t] . 

In Fig. |l^ we show lines for this alternative fit together with measured data. Due to the presence of the /i^^f term, 
fit lines show a small increase close to the chiral limit of valence quark when the difference between sea and valence 
quark is large. This is similar to the behavior observed for quenched QCD in Ref. Q. The amount of increase becomes 
smaller when sea and valence quarks have values closer to each other, and vanishes for full QCD. 



Fit parameters and x /^df are given in Table VIII . x /^df is slightly smaller for the fit with Eq. (49) than the 



one with Eq. ( [fq ) but the difference between the two is not significant. We can therefore not answer the question 
whether a fit with power 3/2 or 2 is preferred. We employ Eq. (^ ) for main results and use Eq. ( ^ ) to estimate 
systematic differences arising from the choice of chiral fit form. 



C. Baryons 

Baryons are made from three valence qua rks a nd hence their masses are expressed in terms of the three /i^'s and 



/igca- In the measurements described in Sec. II B, however, at least two valence quarks are degenerate. We use /i2 to 
stand for the pair of degenerate valence quarks and /ii for the third valence quark. 

For the decuplet baryons masses can be expressed as a function of the average valence quark mass. Hence we define 

Mval = (Mi + M2 + Ai2) , (50) 

and plot decuplet baryon masses as a function of /ivai in Fig- |2[ The behavior of mass data is very similar to the 
one observed for vector meson masses with clearly distinguishable lines of variable slope for partially quenched data 
and stronger curvature for full QCD data. We therefore employ an ansatz of the same structure as for vector mesons 
which takes the form 

mu ^Ksoa; ^vai' '^val ' '^val) ~ + Msca + -B^/ival + C'^ fl^^^ + fi^^y + C^/iscaAival- (51) 

As shown in Fig. O data are fitted well with x^/^df of at most 0.35. Fit parameters and x^/^df are given in 



Table IX 



Octet baryon masses are not simple functions of the average valence quark mass. This can be seen in Fig. |13| where 
we plot masses of S-like octet baryons as a function of /Xvai defined in Eq. (|5^) . The three sets of partially quenched 
data VVV, SVV and SSV lie along different lines. We also see a clear distinction between results for different sea 
quark masses. 

We analyze octet baryon masses by using a formula inspired by chiral perturbation theory pO[ | . In the leading order 
E-like and A-like octet baryon masses are parametrized as a function of quark masses with two constants bo and 
bo. We use these expressions for terms linear in the valence quark mass. For convenience we use a slightly different 
notation; the parameters and D'^ are related to those of Ref. |5^ through = — 2(6d 4- bo) and D'^ = —2bo. 
In order to describe the dependence on the sea quark mass we add linear terms in the sea quark mass, and terms 
quadratic in the sea and valence quark mass to incorporate curvature seen in mass data. The number of additional 
terms introduced by this procedure is limited by the requirement that m\ — niY: when /ii — ^2- This leads to 
expressions for S-like and A-likc baryons of the form 



(«sca; 4aP 4al> 4al) = + ^^A^sca + {F^ ~ D^) fi, + 2F^ ^i2 

+ C^i^l, + C„o,MiM2 + K + Cf) M? + K - C^) M2 

+ (C£ + CfJ Msoa/il + (C£ - CfJ /iscaAi2, (52) 
mA («.ea; 4ai, 4ai) = + Sf^soa + + ^\ Ml + 2 [f^ _ ^) ^2 
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Figure 13 shows masses and fit for S-hke octet baryons. Different hne styles are used for the three types of partially 
quenched data, VVV, SVV and SSV. They do not fall onto each other because of the presence of monomial terms in 
fii in Eq. (|5^). Fit parameters and x^/^df are given in Table 

D. String tension and Sommer scale 

In full QCD, gluonic quantities are still subject to chiral extrapolations through their indirect dependence on sea 
quark masses. We therefore perform such extrapolations on the parameters descr ibing the static quark potential. 
In Fig. H we show ^/a^ and l/ro obtained from the analysis described in Sec. tllQ, as a function of the squared 



pseudoscalar meson mass with valence quarks equal to the sea quark. The sea quark mass dependence of both 
quantities is approximately linear. Therefore we apply fits of the form 

and 

(l^sca) = ^" -^ro^ps (^sca! '^sca; '^sca) ■ (54) 

for extrapolations to the chiral limit, and I/vq in the chiral limit are given in Table 0. 

V. FULL QCD LIGHT HADRON SPECTRUM 
A. Determination of the physical points 



Using the chiral fits of Sect. |V| we determine the physical point of quark masses and the lattice spacing for each 
p. As experimental input we use = 0.1350 GeV and Mp = 0.7684 GeV for the up-down quark sector. For the 
strange quark sector, we compare the two experimental inputs Mk — 0.4977 GeV and — 1.0194 GeV. 

The two flavors of dynamical quarks in our simulation represent up and down quarks which are taken as degenerate. 
Hence we set /ivai — Msoa in Eq. (|8|) and determine the pion mass m-Tr in lattice units by solving the equation 

^MlL (^K) 

Ay + {BY + BY)ml + {CY + CY + CYJmt M, ^ ^ 

for m^. The rho meson mass in lattice units nip is then found by inserting m^r into Eq. (|4^). The error is determined 
with the jack-knife procedure described in Appendix The result of nip is used to set the lattice spacing a by 



identification with the physical value Mp. Lattice spacings obtained in this way are given in Table XI. Inserting 
obtained just above into Eqs. (|T]) and ( |5^ ) with = Hsca. = iti^ the masses of non-strange baryons N and A are 
determined. 

We calculate the strange spectrum in two ways, using either the mass of X or meson as input. As a preparation, 
we determine the hopping parameter of up and down quarks Kud by solving the equation mpg{Kud', Kud, Kud) 



applying the chiral formula Eq. (40) and substituting obtained above. The hopping parameter corresponding 
to the strange point Ks is then fixed by the relation mpg(K„£i; Ks)/to^ — M'^/M'^. In the next step, Ks is 
used to determine the mass of the rjss, a fictitious pseudoscalar meson consisting of two strange quarks, through 
TO^^^ = TOpg(«;„£j; Kg, Ks)- Finally, values of and m^^^ are inserted into Eqs. (^), ( |5l| ) and js^ ) to obtain the rest 
of the spectrum. 

In an alternative determination using the (j) meson mass as input, we first calculate the mass of the rj^s meson by 
using Eq. (E^) and solving the equation 



+ BY ml + BYm^^^ + CY ml + CYmf^^^ + C.^m^m^^^ _ M, 



for m^i^^. Substituting and m^ the spectrum can be calculated as above, except for the K meson, for which first 
Ks is determined from m^^^ — rripg 

{Kud] Ks, Ks ) and then inserted into = mlg{Kud; Kud, 
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We list latti ce sp acings and the hopping parameters K^d and Kg in Table XL Results for the hadron spectrum are 



given in Table KIl. In Fig. |15| hadron masses are plotted as a function of the valence quark mass fi^ai- For this figure 

Tq is used to plot data at different lattice spacings together. Lines, 



obtained from Eqs. (^8|), ( |5l| ) and (|52), correspond to a partially quenched world with sea quarks equal to the physical 
up and down quarks. 

B. Continuum extrapolation 

In Fig. |l^ we show meson masses as functions of the lattice spacing. Baryon masses are plotted in Figs. |l^ and 
p^ . Solid symbols represent our main results at three lattice spacings with a constant physical lattice size. Additional 
masses a.t f3 — 2.2 with a smaller lattice size are depicted with open symbols. 

We find that scaling violations are contained within acceptable limits. The largest scaling violation for mesons is 
observed in the K meson mass (using (j) as input), which changes by 6% between a = 0.22 fm and a = 0.11 fm. The 
largest difference in baryon masses between these two lattice spacings occurs with A for decuplet baryons and with 
S (with K as input) for octet baryons, both amounting to 3%. 

The RG-improved gluon action leads to scaling violation which starts with O(a^). With our quark action, since 
the clover coefficient csw is not tuned exactly at one- loop order, the leading scaling violation is 0{g^a). Here 
is the renormalized coupling constant g^—[fi) |5l) evaluated at a fixed scale /i, which is a constant. Higher order 

perturbative corrections of order g^aloga can be neglected because in our short range of lattice spacings log a is 
almost constant. Accordingly, we attempt continuum extrapolations by applying linear fits to the main data at three 
lattice spacings. We do not include results at = 2.2 because of its smaller lattice size compared to the other runs. 
Lines from linear fits are plotted in Figs. |l^, ^ and [l^. The slopes of the fits are small; parametrizing the dependence 
on the lattice spacing as m = 'Ticont(l — aa), we find, using tuk as input, typical values of a « 0.02-0.04 GeV for 
mesons, a w —0.005 GeV for octet baryons and a « 0.04-0.07 GeV for decuplet baryons. 

The values of for these fits are x^/Nuf ~ 5-7 for mesons, resulting in a goodness of fit Q w 1-2%. The quality of 
fits is therefore marginal. Partly due to larger error bars, fits for baryons are better with /^bf ~ 2 corresponding 
to Q ~ 15%. Having only three data points does not allow us to explore the magnitude of possible higher or der t erms 



of scaling violations. Hadron masses extrapolated to the continuum limit with linear fits are listed in Table XH 



C. Hadron spectrum in the continuum limit 

We observe that meson masses in the continuum limit are quite close to experiment. When using K as input, the 
differences for K* and (j) are 0.7% or 1.3%, respectively, which amount to 1.6 a or 1.9 a in terms of the statistical 
error. When using (j> as input, the mass of the K* is within 0.2% of experiment while the K mass differs by 1.3% 



which is still within the statistical error. As we discuss in more detail in Sec. VII, these results are markedly improved 
from those of quenched QCD ||] which show deviation of about 10% from experiment. 

The situation is different for baryon masses. As is seen with S and E in the octet in Fig. ^ and with in the 
decuplet in Fig. |l^, there is good agreement with experimental masses when the strange quark content is high. The 
difference from experiment increases as strange quarks are replaced with up-down quarks, and the largest difference 
is observed for non-strange baryons; the nucleon mass is larger than experiment by 10% or 2.6 a, and the difference 
for the A is 13% or 2.8 a. 

This pattern of disagreement with experiment appears to be present already at finite lattice spacings. Hence it is 
likely to be a systematic effect rather than a statistical fiuctuation. A possible reason behind this are finite size effects 
arising from the lattice size of La ~ 2.5 fm. We expect lighter baryons made of lighter quarks to be affected more 
from these effects, which is consistent with the pattern we observe. A detailed investigation is needed, however, since 
finite size effects in full QCD can be quite complicated, arising from both sea and valence quarks wrapping around 
the lattice in the spatial directions. 

We add a remark for strange baryons. Masses obtained using either K or as input (left and right panels in 



Figs. [17| and |18|) differ at coarse lattice spacings. The difference decreases with lattice spacing, however, and almost 
disappears toward the continuum limit. This reassuring finding is connected with a good agreement of the strange 
meson spectrum with experiment in the continuum limit. 

VI. QUENCHED QCD WITH IMPROVED ACTIONS 
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A. Purpose 



Up to this stage we have discussed the two-fiavor full QCD hadron spectrum. In order to analyze how dynamical 
sea quarks manifest their presence in the spectrum, we need to compare full QCD results with those of quenched 
QCD. 

The quenched hadron spectrum has been examined in detail in Ref. Systematics of simulations in Ref. |^ differ, 
however, from those of two-flavor QCD in the present work. The standard plaquette gluon action and the Wilson 
quark action are used in Ref. and the continuum extrapolation is made from a finer range of lattice spacing 
a « 0.1-0.05 fm in ||^ as compared io a ~ 0.2-0.1 fm in the present work. The lightest valence quark mass is pushed 
down to TOps/my ~ 0.4 for quenched QCD while it only reaches mpg/my « 0.5 in full QCD, and the physical lattice 
sizes are La « 3 fm for quenched QCD and La « 2.5 fm for full QCD. 

We consider that a more direct comparison with a common choice of actions over similar range of lattice parameters 
is desirable. Therefore we carry out a new set of quenched simulations with the same set of improved actions as 
employed for two-flavor full QCD. 



B. Matching quenched and full QCD simulations 

We use the string tension to match the scale of quenched QCD with that of full QCD, i.e., for each value of (3 and 
Ksoa at which full QCD simulations are made, we make a corresponding quenched run with f3 chosen such that the 
string tension cr in lattice units takes the same value. 

This is carried out at four values of Ksea at P = 1.95 and at 2.1, and also at the chiral limit Kgca = i^c at the two 



values of 13 of full QCD. A summary of the 10 gauge couplings used for quenched simulations is given in Table XIII. 
In the same table we list measured string tensions, to be compared to the ones for full QCD in Table We also 
quote lattice spacings obtained using the rho meson mass as input. 

Simulations are carried out using the same lattice size as the corresponding full QCD runs, namely 16"^ x 32 and 
24^ X 48. Physical lattice sizes vary therefore between La ~ 2.6 fm and La « 3.5 fm. 



C. Simulation details 



Gauge configurations are generated with a combination of the 5-hit pseudo-heat-bath algorithm with two SU(2) 
sub-matrices and the over-relaxation algorithm. The two algorithms are mixed in the ratio of 1:4 and the combination 
is called an iteration. For vectorization and parallelization of the simulation code, a 16-color algorithm is developed 
for the RG-improved gauge action. 

We skip 100 iterations between two configurations for hadron propagator measurements. We check that this number 
of iterations is sufficient to regard the configurations to be independent. We calculate hadron propagators over 200 
configurations per gauge coupling. These statistics are comparable to the number of independent configurations in 
the fuU QCD runs. 

The measurement procedure parallels the one for full QCD. Hopping parameters are chosen so that ratios mpg/my 
for degenerate mesons match the ones of the corresponding full QCD run. For the quark matrix inversion we use the 
same set of stopping conditions and smearing parameters as the ones for corresponding full QCD runs. Masses are 
extracted from hadron propagators with smeared sources using correlated fits and fit ranges similar to those used for 
full QCD. 

For chiral extrapolations we follow the strategy of fitting vector and baryon masses as a function of measured 
pseudoscalar masses, and these in turn as a function of valence quark masses. To be specific, we fit pseudoscalar 
meson masses to the formula 

where variables are defined as in Eqs. (^8|) and (^9|). This is the quenched analogy of Eq. (|4^) with terms containing 
"^sca'^ dropped. Similarly, when making fits as a function of AWI quark masses we employ the formula 

™ls 4ai) = Km^T + c: {m^TY , (58) 

which corresponds to Eq. (^) for full QCD. 
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For vector mesons an inspection of mass data, plotted in Fig. shows that they are weh described by a Hnear 
function. If we nevertheless perform a quadratic fit the coefficient of the quadratic term is ill defined with large error 
bars. We therefore employ fits with 

"^V , = + ^Mval , (59) 



as shown in Fig. [19|. Parameters of chiral fits to mesons with Eqs. (^Tj), (^8|) and ( p39D are given in Table XIV 



Analysis of baryon masses proceeds in a similar way. For decuplet baryons we again find quadratic terms in quark 



masses to be unnecessary. Data for baryon masses together with fits are plotted in Figs. 20 and ElJ 



D. Results 



Physical hadron masses are summarized in Table XV. They are plotted as a function of the lattice spacing in 
Fig. for mesons and in Figs. ^ and for baryons. 

In the same figures we also plot hadron masses obtained with the standard action in Rcf. 1^. In this work, the 
analysis was made with two sets of functions for chiral extrapolation. The main method used functional forms 
predicted from quenched chiral perturbation theory. As an alternative method polynomial fits were also employed. 
It was found that results from the two methods are consistent with each other within errors after the continuum 
extrapolation. In particular, conclusions on the deviation of the quenched spectrum were not altered by two different 
methods. Since in this work we use polynomial fits for the analysis, we take hadron masses from polynomial fits in 
Ref. d for a comparison within quenched QCD. 

We perform continuum extrapolations of hadron masses fo r the improved action linearly in the lattice spacing in 
accordance with the leading scaling violation discussed in Sec. |VB| . Good x^/^df ~ 1 are obtained for meson masses. 
Baryon mass data exhibit some scatter and as a result larger x'^Z-^df are observed. The largest value, reached for 
the E baryon, is x^/-^df — 2.8; hence we consider the scatter to be still within the limits of statistical fluctuations. 

Comparing masses in the continuum limit a good agreement is found between calculations with the standard and 
improved actions. All results are consistent within the statistical accuracy. This is a confirmation that the quenched 
light hadron spectrum deviates from experiment |^ . 

Meson masses from the two choices of actions both show very good scaling, and they are already in agreement even 
at finite lattice spacings. For baryons scaling behavior is improved for the improved action. This is in accordance 
with our initial study of action improvement [ p2| , notwithstanding that this study was carried out for full QCD. The 
largest scaling violation in improved baryon masses is observed for the nucleon with a difference of 14% between 
~ 1 GeV and the continuum limit. 

VII. SEA QUARK EFFECTS IN THE LIGHT HADRON SPECTRUM 

A. Light meson spectrum 



In Fig. 25 we compare the continuum extrapolation of vector meson masses using the K or cj) meson mass as input 
for full QCD and for the two quenched calculations. The deviation of the quenched spectrum from experiment is 
considerably reduced in full QCD. For the K* meson the deviation is reduced from 2.6% (3.1% with the standard 
action) to 0.7%, and for the meson from 4.1% (4.9%) to 1.3%, if the K meson mass is used as input. Using the 
(j) meson mass as input, the difference in the K* meson is less than 1% for both quenched and full QCD, while the 
deviation for the K meson is reduced from 8.5% (9.7%) in quenched QCD to 1.3% in full QCD. We consider this 
improvement in the meson spectrum to be a manifestation of sea quark effects. 

An important factor in reaching this conclusion is the continuum extrapolation. At finite lattice spacings the 
difference between full and quenched QCD is not obvious. At two coarse lattice spacings in particular, the two sets 
of data are roughly consistent. However, the trend towards the continuum limit is different. Full QCD leads to an 
increase for the K* and meson mass (decreasing for the K meson mass) in contrast to a flatter behavior in the 
quenched masses. A support that these trends are not just fluctuations is provided by the additional calculation at 
f3 = 2.2, showing higher (lower) lying values, as can be seen from small flUed circles in Fig. 

Let us discuss systematic errors which are relevant for this conclusion. In Fig. ^ we show how the K* meson mass 
changes when different functional forms are used for chiral extrapolation. Filled squares represent masses obtained 
using the flt with Eq. ( p9| ) instead of our standard analysis plotted with filled circles. There is a noticeable effect on 
the K* mass, which increases by 1% in the continuum limit. A similar effect is seen for the quenched data where we 
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show results of Ref. |5[] for two ways of chiral extrapolation. The trend remains, however, that the continuum value 
for full QCD lies much closer to experiment than in quenched QCD. 

Another source of systematic errors is the continuum extrapolation. Within the small number of data points 
available for full QCD, we may estimate the upper error by making an extrapolation from the two points at [3 — 1.95 
and 2.1, and the lower error by taking the value aX [3 — 2.1. For the K* meson mass this yields rriK* — 0.890(4)^9''^ GeV 
where the second error represents the systematic error estimated in this way. For a complementary analysis in the 
quenched simulation with the improved action, we make a linear fit to the five points with fine lattice spacings 
corresponding to the full QCD point at /3 = 2.1 for the upper error, and take the left-most point with the finest lattice 
spacing for the lower error. We then obtain ttik* = 0.873(2)12 GeV. 

Similar analyses lead to = 1.007(7)^^7 GeV and ttik — 0.504(8)^25 GeV for full QCD compared to = 
0.977(4)1^4^ GeV and rriK = 0.540(5)j^^g GeV for quenched QCD with improved actions. Hence systematics of the 
continuum extrapolation are unlikely to annul a closer agreement of full QCD masses with experiment compared to 
quenched QCD. 

In summary we find that effects of dynamical sea quarks are present beyond the systematic as well as statistical 
uncertainties in strange meson masses. 



B. J parameter 

A useful quantity to quantify sea quark effects in the meson sector is the J parameter ||5^ defined by 



(60) 

mv/mps=MK* /Mk = 1.S 



where only valence quark masses are to be varied in the differentiation. In the real world this corresponds to a 
comparison between strange and non-strange mesons. The derivative in Eq. ( |60| ) can be replaced by a finite difference 
and an "experimental" value for J is then obtained as 

We calculate J from fits to vector mesons as functions of pseudoscalar mesons in two different ways. In the first one 



we use combined fits with Eq. (48), keep /igca fixed and calculate derivatives with respect to /ivai- This leads to the 
curves shown on the left side of Fig. |2^. For the second method we employ separate partially quenched fits for each 
simulated sea quark. We use quadratic fit functions obtained from dropping all terms containing /igoa in Eq. (|4^). 
Results are plotted with filled symbols in Fig. |2^. They tend to scatter more since, in contrast to combined fits, no 
smoothness in the sea quark mass is imposed for separate fits. The two methods yield consistent results within at 
most two standard deviations, showing a trend of increase as the lattice spacing is reduced. At fixed lattice spacing, 
on the other hand, we do not see a clear dependence as function of the sea quark mass. 

On the right hand side of Fig. ^ we plot J at the physical point for quenched and two-flavor full QCD as a function 
of lattice spacing. For quenched QCD, the values do not show much variation, and a linear extrapolation to the 
continuum lim it give s J = 0.375(9)j]2^ where the second error represents the systematic error estimated in the same 



way as in Sec. VII A, This is consistent with earlier observations of a too small value of J in quenched QCD. 

Full QCD data at /3 = 1.8 and 1.95 do not differ much from this value. It is intriguing, however, that at /3 = 2.1 
(and also /3 = 2.2) J is sizably larger. Consequently the continumn value of J = 0.440(15)^27, estimated by a linear 
extrapolation, lies much closer to experiment. 



C. Sea quark mass dependence 

An interesting question with dynamical sea quark effects is how their magnitude depends on sea quark mass. We 
examine this point by calculating the mass ratio rriK' / "m-p for fixed valence quark masses as a function of sea quark 
mass. 

The analysis proceeds in the following steps. We leave the sea quark mass parametrized by ^sca as a free parameter, 
and first determine the valence pion mass "tott" and the rho meson mass "mp" corresponding to a given ratio 
TOps/mv ="''7T-7r"/"™p" which may be different from the physical one, e.g., tops/tov = 0.5 in an example shown below. 
In the next step the strange pseudoscalar meson mass "to^,," is fixed by a phenomenological ratio "TO,,^^"/"m0" = 
^y2M'^ — M^/M^ — 0.674. To be specific, for full QCD an interpolation to this ratio consists of solving the equation 
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0.674 (62) 



for "m,j^J'. Finally using "m^" and "m,;^/' determined above, and setting /ivai = ("mj" + "to^^^")/2 in Eqs. (|4q ) 
or ( [59I ) we obtain the mass "mjc*" of a fictitious if* meson. In this setup "mp" is again used to set the scale 
by calculating the mass ratio "m^. " / "mp" . As a measure for the lattice spacing "rrip" in lattice units is used for 
continuum extrapolation. 

In Fig. ^we illustrate the ratio "m;f"/"TOp" as a function of mps/mv of sea quarks when rnpg/mv of the valence 
quarks is fixed to 0.5. Naively we would expect the points to be a smoothly decreasing function of mps/mv, reaching 
the quenched value at mpg /my = 1 corresponding to infinitely heavy sea quark. In contrast to this expectation, but 
consistent with the findings for the J parameter, sea quark effects are almost constant up to mpg/my « 0.7-0.8, 
which roughly corresponds to the strange quark. This may be an indication that sea quark effects turn on rather 
rapidly when sea quark mass decreases below a typical QCD scale of a few hundred MeV. 

VIII. LIGHT QUARK MASSES 

Hadron mass calculations in lattice QCD provide us with a unique and model-independent way to obtain quark 
masses. The main findings of our light quark mass calculation have been presented in Ref. |3^]. We give here a more 
detailed account of the analysis and results. 

A. Extraction of quark masses 

Quark masses can be calculated by inverting the relation (Eo) and (p4|) between quark masses and pseudoscalar 



meson masses, and substituting and m^^ determined in Sec. V A 

.(1) _ ^(2) 
"val ~ '^val 



For the average up and down quark mass, we set k^J, = n'^} = ^sca and evaluate the hopping parameter K^d 



vwi _ (l/Kud — l/'*c)/2 where Kc is the critical hopping parameter where the pseudoscalar meson mass made of 



for these quarks by solving the equation TOpg(K 

ud: f^ud: ^ud) — The V^VI quark mass is then determined by 

^ud 

sea quarks vanishes mps(Kvai = ^sca. = i^c) = 

An alternative definition for the VWI quark mass, called partially quenched VWI quark mass (VWI,PQ), has been 
proposed in Ref. |5^. The partially quenched (PQ) chiral limit is defined as the point of Kvai where the pseudoscalar 
meson mass vanishes for fixed Ksea, and the corresponding hopping parameter is denoted as t^'^- As apparent from 
Fig. |9[ values of k^*^ exhibit a clear dependence on Kgoa and coincide with Kc only in the limit Kgea = Kc- The 
proposal in Ref. |Q consists of defining the quark mass via m^^"^^ = {l/nud — ^/k^'^)/2 where for the value at 
Ksca = Kud is substituted. This is equivalent to a fictitious situation where the simulation is performed with dynamical 
quarks at their physical value of up and down quarks, the spectrum of pseudoscalar mesons is measured for several 
values of the valence quark and the chiral limit is defined at the point where masses of pseudoscalar mesons vanish. 

A third determination of the average up and down quark mass is obtained using the AWI definition of quark mass. 
It is unambiguously determined from Eq. ( ^ ) by setting m^^^ = m^^^ and solving for m^g — m^. 

The determination of the strange quark mass is made in a similar way. Keeping the sea quark mass fixed at the 
average up and down quark mass determined above, i.e., Ksca — n^d in Eq. (^) and m^^^ = m^^ in Eq. (^4|), we 
calculate the point of strange quark by tuning Kvai or m^^^ so that Tipg equals m^^ obtained from the spectrum 
analysis. 

Since m^^^ depends on the physical input, the strange quark mass also depends on this input, and we consider the 
two cases where the K meson mass and the (j) meson mass are used as input. In an exact parallel with the average 
up and down quark mass, we calculate the strange quark mass with three definitions. 

Bare quark masses are converted to renormalized quark masses in the modified minimal subtraction (MS) scheme at 
/i = 1/a by the use of one-loop renormalization constants and improvement coefficients, summarized in Appendix |^. 
For the two definitions of VWI quark mass this consists of a conversion of the form 

/ mVWIX ™VWI 

mvwi ^z^(l + b^"^ ) ^ , (63) 

while the renormalized AWI quark mass is obtained with 
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ignored this contribution. After conversion to the MS scheme we employ the three-loop beta function to run quark 



Since (hA — hp)m^^^ ju^ ~ — 0.0019(j(|jgTO^^V'"o ~ 0.0006 ^ 1 is negligible even for the strange quark, we have 



masses to a common scale ol [i — 2 GeV. Numerical results are listed in Table XVI 



B. Continuum results and systematic uncertainties 

Quark masses are plotted as a function of the lattice spacing in Figs. ^ and In these figures we also show lines 
for continuum extrapolations performed for each definition of quark mass separately. For extrapolation we employ fits 
linear in the lattice spacing, corresponding to the leading order scaling violation. We only include data from runs at 
three lattice spaci ngs fo r extrapolation, leaving the run at /3 = 2.2 for a cross-check. Results of these extrapolations 



are given in Table XVI 



For rriud scaling violations are very small if the AWI definition is used. The difference between the value at the 
coarsest lattice spacing and the continuum value from a linear extrapolation is only 1%. In contrast, the two other 
definitions show sizable scaling violations. The partially quenched quark mass at the coarsest lattice spacing is 20% 
higher than in the continuum limit while the VWI quark mass is lower by 34%. Furthermore, the VWI quark masses 
exhibit some curvature. 

The situation is similar for the strange quark mass when the K meson mass is used as input. Scaling violations are 
small for the AWI quark mass, amounting to a value 3% higher at the coarsest lattice spacing than in the continuum 
limit. For the two VWI quark masses, on the other hand, this difference amounts to 15%. If the meson mass is 
used as input scaling violations are larger. In this case even the AWI quark mass is 30% larger at the coarsest lattice 
spacing than in the continuum limit and for the two VWI quark masses the difference is as large as 45%. 

Having data at only three lattice spacings, it is difficult to explore scaling violation for each definition of quark 
mass in further detail. An important observation for linear continuum extrapolation is the fact that the different 



fits to each definition converge in the continuum limit within two-sigma of statistics (see Table XVI). In particular, 
VWI quark masses, where the largest scaling violations are observed, are consistent with AWI masses, where scaling 
violations are generally small. This leads us to perform a further fit, linear in the lattice spacing and having a common 
continuum value. With such fits we obtain mud = 3.44(9) with x^/^df = 2.9 and nis = 88.3(2.1) with x^/A^df = 1-3 
{K input) or = 89.5(4.3) with x^/-^df = 3.0 {4> input). These masses lie between the ones from individual fits 
and can be considered as a weighted average. We utilize these numbers for central values of quark masses. 

The errors quoted above are only statistical. Systematic errors arise from the continuum extrapolation, the chiral 
extrapolation at each lattice spacing, and from the use of one-loop renormalization factors in relating the lattice values 
of quark masses to those for the continuum. 

One way to examine systematic errors in the continuum extrapolation is to include higher order terms in the 
combined fits. Such fits, however, are unstable and do not lead to higher confidence levels, in particular for mud- We 
therefore estimate uncertainties of the continuum extrapolation from the spread of values obtained by separate fits to 
data from the three definitions. Taking differences between the values from separate fits and that from the combined 



fit leads to the errors quoted in Table XVII 



We estimate the error from chiral extrapolation by changing the fit formula. The functional form used for the 
determination of physical points, and hence quark masses, is given with Eq. (^8|). Changing this to the alternative 
form of Eq. ( ^ ) has several effects which, combined together, lead to a decrease of the continuum value by 2-8% 
from the main analysis. This is used as estimate of the lower error. For the upper error we add cubic terms m^ to 
the formulae (^0|) and (^J) for pseudoscalar mesons as functions of quark masses. This results into an increase of the 
quark masses at each lattice spacing and also in the continuum limit. 



Turning to the problem of renormalization factors, we list one- loop corrections in Table XXVI. Their contribution 
is at most 13% at the strongest coupling, and hence we may expect higher loop contributions to be smaller. Since a 
non-perturbative determination of the renormalization factors is yet to be made for our improved actions, we estimate 
effects of higher order corrections with a shift of the matching scale from n = 1/ a to fj, = tt/o, and with use of an 
alternative definition of the coupling given in Eq. (p3[). The former leads to an increase by 2%, while the latter leads 
to a decrease of 5-7%. 



Finally we add the statistical and the systematic errors listed in Table XVII in quadrature to obtain the total error. 
This leads to the final values 

m^(2GeV) = 3.44ll];i4 MeV, (65) 
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for the average up and down quark mass and 



mf^(2GeV) = 88^^ MeV input, (66) 

= 90tli MeV input. (67) 



for the strange quark mass. 



C. Sea quark effects on light quark masses 

In Figs. O and we compare quark masses in full QCD (filled symbols) with those in quenched QCD (open 



symbols). The quenched results for improved actions (thick open symbols) are obtained from the analysis of Sec. VI 
in parallel to those of full QCD. There is no ambiguity in choice of the critical hopping parameter and so there is only 
one definition of VWI quark mass. We also show quark masses for the standard action reported in Ref. |^ (thin open 
symbols). 

Long dashed lines are from the combined fit for full QCD, for which the errors drawn in the continuum limit 
includes the systematic errors. The continuum limits for quenched QCD are estimated with a combined linear 



continuum extrapolation. They are listed together with quark masses in full QCD in Table XVIII. 

Comparing the two quenched calculations of quark masses we first note that scaling violations are visibly reduced 
for the improved action. This is most noticeable for the strange quark mass where masses from improved actions show 
a flat dependence against the lattice spacing a, while they exhibit a sizable slope for the standard action. Nonetheless 
quark masses in the continuum limit from the two calculations are in good agreement. 

This confirms an inconsistency of 20-30% in the quenched estimate of the strange quark mass depending on 
whether the K meson mass or the (j) meson mass is used as input. 

A comparison of full and quenched QCD establishes that the effect of dynamical quarks decreases estimates of 
quark masses. This point was previously argued from renormalization-group running of the gauge coupling and quark 
masses in Ref. pst . For two dynamical flavors examined in the present work mud becomes smaller by about 25%. For 
the strange quark the decrease is 20-25% using K as input, and 30-35% for cj) as input. 

In two-flavor full QCD the strange quark mass is consistent between the two different inputs within the errors of 
5-10%. This is caused by a different amount of decrease between quenched and full QCD. Thus the inconsistency 
in the strange quark mass of quenche d QCD almost disappears in the presence of two flavors of sea quarks. This is 



directly related to the finding in Sec. |VII A| that the K — K* and the K — (p mass splittings show a close agreement 



with experiment while there is a clear discrepancy for quenched QCD. 

IX. DECAY CONSTANTS 

A. Pseudoscalar meson decay constants 

The pseudoscalar decay constant /ps is defined from matrix elements of the axial vector current through the relation 

(OIA4IPS) = /psmps- (68) 

We include the 0{a) improvement term in the axial vector current, and employ one-loop renormalization constants 
as described in Appendix The decay constant is evaluated from the formula 



/ps = 2^uoZa (1 + bA-) ^J^. (69) 

V UqJ C%\ TOPS 

Here for m we substitute the VWI,PQ quark mass, superscripts I and s distinguish local and smeared operators, 
and various amplitudes are extracted in the following steps. The pseudoscalar mass tops and the amplitude Cp are 
determined from 

(P'(t)P^(O)) = [exp(-TOpsi) + exp(-mps(Lt - t))] . (70) 
Values of tops are listed in Appendix 0. Keeping the mass fixed, we extract Cp and C% from the fits 
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(P'(i)P'(O)) 
{A[{t)P%0)) 



&P [exp(- 
C% [exp(- 



-mps<) 



exp(- 
exp(- 



mps{Lt 
rnpsiLt 



t)]. 



(71) 
(72) 



The chiral extrapolation of the decay constant is carried out in the same way as for vector meson masses. Hence 
we employ a combined fit in sea and valence quarks of the form, 



f f (1) 

/pS yKsca,; k'J^, K. 



2 

sea 



+ /iyai + C^/iscaMval, 



(73) 



where the /i's have the same meaning as in Sec. EVB 
are shown in Fig. 



Sec. VA 



Parameters of the fit are given in Tabic |XIX 
) Msoa = ml and /ivai = ml or /Xvai 



Pseudoscalar decay constants together with fits with Eq. ( [T^ ) 



Setting in Eq. (rZilJ /Zgoa = ™^ and /ivai = rn 
/tt and /if are obtained in lattice units 



i^)/2 obtained from the spectrum analysis in 



Decay constants in physical units are finally calculated using the 



XX 



lattice spacing from the rho meson mass and are listed in Table 

The extraction of decay constants in quenched simulations is made similarly to that in full QCD. For chiral ex- 
trapolation a simpler version of Eq. (^3|) ignoring sea quark mass dependence is used, and the quadratic term /i^^^j 
is dropped as linear fits already yi eld g ood as illustrated in Fig. /^r and fx obtained from calculations in 
quenched QCD are quoted in Table XX . 

In Fig. ^ we show the lattice spacing dependence of and Jk in full and quenched QCD. For a comparison 
we also include results obtained in quenched QCD with the standard action ||^. The most noticeable point is large 
violation of scaling in full QCD. The values at the coarsest lattice spacing a = 0.22 fm are 50% larger than that at 
the finest lattice spacing of a = 0.11 fm. Scaling violation is milder for quenched QCD, but still decay constants at 
a — 0.22 fm are 15% larger than those at a = 0.11 fm. 

The origin of large scaling violation in the pseudoscalar decay constant is not clear at present. Possible origins are 
contributions of higher order corrections in the renormalization factors and 0(a) terms in the axial vector currents. 
A suggestive hint pouting toward these origins is provided by the ratio Jk/ fir — 1 for which such corrections may 
largely cancel out. As shown in Fig. |3^, one observes much reduced scaling violation for this quantity. Furthermore, 
a trend of increase toward the experimental value as effects of sea quarks are included is also apparent. 



B. Vector meson decay constants 



Vector meson decay constants are defined as 



(0|F,|V) = e.Fvmv, 



(74) 



where ti is a polarization vector and my is the mass of the vector meson. 

The numerical procedure employed to calculate vector meson decay constants parallels the one for pseudoscalar 
decay constants. As discussed in Sec. [II A, the rho correlator with smeared source is fit with 



{V\t)V''{Q)) = C'y [exp(-TOvi) + exp(-mv(it - t))] 
which determines my and Cy. Using my as input we make fits to the correlator 

{V\t)V\{))) = C'y [exp(-myt) + exp(-mv(it - t))] , 



(75) 



(76) 



where the amplitude Cy is the only fit parameter. Renormalized vector meson decay constants are then obtained 
through 



(77) 



Fy = 2kwo^v 1 + hy 




where expressions for pcrturbative renormalization factors are given in Appendix and for m we su bstitu te the 
VWI,PQ quark mass. We note in passing that we do not include the improvement term cyduTn^ii, in Eq. ( C15| ), since 
the corresponding correlator has not been measured. 

For chiral extrapolations we again employ combined quadratic fits a s defined by Eq. (^3|) . These fits describe the 
data well, as shown in Fig. Fit parameters are given in Table XIX. Vector meson decay constants obtained from 
quenched simulations are plotted in Fig. |3^. As for pseudoscalar decay constants they are well described by linear 
fits. Final values of Fp, Fk* and F^ in physical units are listed in Table XX for both full and quenched QCD. 
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The lattice spacing dependence of Fp and in full and quenched QCD is shown in Fig. |39|. We again include results 
obtained in quenched QCD with the standard action [|| for comparison. Vector meson decay constants in full QCD 
exhibit scaling violations similar to those found for pseudoscalar decay constants; e.g., Fp is 40% larger at a = 0.22 fm 
than at a = 0.11 fm. Consequently, a continuum extrapolation poses similar difficulties as for pseudoscalar decay 
constants. 

Since scaling violation is similar in vector and pseudoscalar decay constants, one may examine the ratio of Fp 
to /tt. The lattice spacing dependence is much reduced for this quantity (see Fig. |40| ), and i^p/J^r is consistent with 
experiment within the error of 5-10%. In contrast to pseudoscalar decay constants, sea quark effects are not apparent. 



C. Non-perturbative renormalization factors for vector currents 

For the clover quark action one can define a conserved vector current which reads 

VFin) - i {^n+^.Ul^{l^ + l)9n + 7 n^uA^^ " l)5n+/,} • (78) 

The non-renormalization of this current can be used to obtain a non-perturbative estimate of the renormalization 
constant for the local current mJsS] according to the relation. 



The non-perturbative renormalization factors obtained from Eq. ( [79| ) and extrapolated to zero quark mass are 
plotted as a function of the gauge coupling constant in Fig. B^. In the same figure we also plot mean-field improved one- 
loop perturbative renormalization factors as calculated in Appendix Non-perturbative values of Zy are significantly 
smaller than those obtained from perturbation theory. This may be partly due to corrections of 0{a) which are 
necessarily included in Zy calculated from Eq. ( [79| ) p^-p6| . 

In Fig. we compare Fp determined with either perturbative or non-perturbative renormalization factors. We 
observe that decay constants calculated with Zy^ exhibit a much flatter behavior as a function of the lattice spacing. 
We take this as an encouraging indication that a further study with non-perturbative renormalization factors will 
help moderate an apparently large scaling violation in the pseudoscalar and vector decay constants. 



X. CONCLUSIONS 



In this article we have presented a simulation of lattice QCD fully incorporating the dynamical effects of up and 
down quarks. A salient feature of our work, going beyond previous two-flavor dynamical simulations, is an attempt 
toward continuum extrapolation through generation of data at three values of lattice spacings within a single set of 
simulations. In order to deal with the large computational requirement that ensue in such an attempt, we have used 
improved quark and gluon actions. This has allowed us to work with lattice spacings in the range a « 0.22 — 0.11 fm, 
which is twice as coarse as the range a<0.1 fm needed for the standard plaquette gluon and Wilson quark actions. 
Still, this work would have been difhcult without the CP-PACS computer with a peak speed of 614 GFLOPS. With 
a typical sustained efficiency for configuration generation of 30-40%, the total CPU time spent for the present work 
equals 415 days of saturated use of the CP-PACS, of which 318 days were for configuration generation and 84 days 
for measurements. 

A major physics issue we addressed with our simulation was the origin of a systematic discrepancy of the quenched 
spectrum from experiment Our new quenched simulation employing the same improved actions as for full QCD 
has quantitatively confirmed the results of Rcf. for both mesons and baryons. 

For mesons, masses in two-flavor full QCD become much closer to experiment than those in quenched QCD. Using 
the K meson mass to flx the strange quark mass, the difference between quenched QCD and experiment of 2.6lQg% 
for the K* meson mass and of 4.ll;" g% for the meson mass is reduced to 0.7^\'l% and 1.312 5% ^^^^ QCD. 
When the (p meson mass is used as input, the difference in the K* meson mass is less than 1% for both quenched and 
full QCD, while the deviation from experiment for the K meson mass is reduced from 8.5j^3 g% in quenched QCD to 
in full QCD. Similarly the J parameter takes a value J = 0.440;^ |^^i in two-flavor full QCD, which is much 
closer to the experimental value J « 0.48 compared to J = 0.375lQQQg in quenched QCD. We take these results as 
evidence of sea quark effects in the meson spectrum. 
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A common point in reaching this conclusion is the importance of continuum extrapolation. Differences between 
quenched and full QCD meson masses are less obvious at finite lattice spacings but the slope of the continuum 
extrapolation is different between them. Unexpectedly, the scaling violation for full QCD is apparently larger than 
for quenched QCD with the same improved actions. A possible origin of this feature is the common choice of csw we 
made for the two cases while the correct csw necessary to remove 0{a) scaling violations need not be the same. 

Full QCD baryon masses exhibit the pattern that the difference from experiment increases with decreasing strange 
quark content. While masses of S and Q are in agreement with experiment, the nucleon mass differs most from 
experiment among the octet, being larger by 10%, and the A among the decuplet by 13%. This pattern of disagreement 
suggests that finite-size effects sizably distort light baryon masses for an La sa 2.5 fm spatial size employed in our 
study. We leave detailed finite-size analyses in full QCD for future investigations, however. 

The sea quark effects in the meson sector have an interesting consequence that the light quark masses decrease by 
about 25% in two- flavor full QCD compared to quenched QCD. An inconsistency of 20-30% in the strange quark 
mass for quenched QCD, depending on the particle used as input, disappears in full QCD within the errors of 5-10%. 

In contrast to the encouraging results above, meson decay constants exhibit large scaling violations which obstruct 
a continuum extrapolation. We have found this trend to be common through light pseudoscalar and vector decay 
constants of this work as well as in heavy- light decay constants . Possibly this problem arises from two- loop and 

higher order corrections in the renormalization factors not included in our analyses. An indication for this explanation 
is given by a much flatter behavior of vector meson decay constants when using a non-perturbative renormalization 
factor derived from a conserved vector current. 

While we consider that the present work has brought sizable progress in our effort toward fully realistic simulations 
of QCD, it is also clear that a number of gaps have to be filled in future studies. One of them is an examination 
of finite-size effects, particularly for baryons. Another is the exploration of lighter values of sea quark masses below 
wps/wv ~ 0.6 for better control of the chiral extrapolation, and generation of data at more points in the lattice 
spacing for a better control of continuum extrapolations. Important in the latter context will be the use of non- 
perturbative improvement coefficients and renormalization factors. Finally, the inclusion of a dynamical strange 
quark will be necessary to remove the last uncontrolled approximation. 
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APPENDIX A: HADRON MASSES 



In Tables XXI - XXIIl we set out the hadron masses measured in full QCD simulations. We list fitting ranges, 

/2 / ]\T-^ -r^ o n /-I T-v-i 00 0(^0 m lo-H-nr»(3 11-1-114-0 fV^iT- oil ttoIiij^iC' r^-f ntn r\ oil r>j-\i-vTl-\i no+-n^T-iC' r-\f i-- o n r\ 1^^^^ C^UOtiG GITOrS 



X /-/Vdf and masses in lattice units for all values of /? and all combinations of k. 
determined with the jack-knife method with a bin size of 10 configurations or 50 HMC trajectories. 



sea and K^^j . 



APPENDIX B: JACK-KNIFE ANALYSIS FOR FULL QCD SIMULATIONS 

In quenched simulations masses of hadrons with different quark content are obtained from the same gauge config- 
urations and are therefore correlated. Often the quality of data does not allow a correlated chiral extrapolation and 
it is usual practice to resort to uncorrelated fits. By using the jack-knife method errors of fit parameters can still be 
correctly determined. 

At first sight the situation seems simpler for full QCD with valence quarks equal to sea quarks. Separate runs 
have to be made for different sea quark masses, and are manifestly uncorrelated. Errors on parameters of chiral 
fits can be correctly calculated from an uncorrelated fit. Nevertheless, the jack-knife method is extremely useful 
even in this case. Since the fit parameters are often highly correlated, the determination of the error of derived 
quantities can not be made with naive error propagation. The jack-knife method takes such correlations into account 
correctly. Moreover, in the setup of two-fiavor QCD entire sets of hadron masses with different valence quark content 
are measured on the same configurations created with a given sea hopping parameter. Combined fits according to the 
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method of Sec. IV have correlations between some of the data, and therefore one is in a similar situation to quenched 
QCD. 

A difference from quenched QCD is that there are as many sets of gauge configurations as sea quarks in the simula- 
tion. They arc mutually independent and can differ in number between runs with various sea quarks. A generalization 
is implemented in the following way. First, hadron masses are determined with the usual jack-knife method. This 
yields mass estimates 'm-\^'^^ (k^^^) for each jack-knife ensemble obtained by omitting the gauge configuration number 
i from the run with sea hopping parameter number k. Mean values and variances are defined by 

-, Nk 

"^H(0 = ^E™H^''(^-a), (Bl) 
i—1 

(AmH(^L))' = E (™H ^^(O - ™h(acL))' • (B2) 

i—l 

Chiral fits are then carried out by replacing mean values fn}i{K^^^) with jack-knife estimates ^j/'^l^^sea) ^^^r the sea 
hopping parameter number k while keeping masses at all other sea hopping parameters at their mean value. This 
procedure gives error estimates (AP)^ as above where P stands for a fit parameter or a quantity derived from fit 
parameters. Since runs at different sea quarks are uncorrelated, the total error AP is determined by quadratic addition 
(AP)2 = I]j.((AP)fe)2. Errors quoted throughout this paper are determined with this method. 



APPENDIX C: RENORMALIZATION FACTORS AND IMPROVEMENT COEFFICIENTS 

In this appendix we summarize renormalization factors and improvement coefficients used in the calculation of 
matrix elements and quark masses. Perturbative calculations to one loop have been carried out in Rcfs. p5| , ^ . 

For the coupling constant we adopt a mean-field improved value in the modified minimal subtraction (MS) 
scheme obtained in the following way. We start with the one-loop perturbative relation between the bare and MS 
couplings for the RG improved gauge action and the 0(a)-improved Wilson quark action [p7[. 



1 1 ll-i^^f 
^ — +0.1000 + 0mU9Nf + ? ^ 



— + 0.1000 + 0.031497V/ + — I ' log(A<a). (CI) 



The formula is reorganized so that 1 /g^ becomes the coefficient in front of in the continuum limit after the mean 
field approximation. Using the one- loop expressions P = 1 — 0.1402(7^ and P = 1 — 0.2689(?2 for the expectation 
value of the plaquette P = (W^^^) and the 1x2 rectangle R = {W^^'^), we obtain the relation 

1 _coP + 8ciP 0.1006 + o.03149iV/ -I- ^^~|^^ og(Ma). (C2) 



Tadpole-improvement is realized by using nonperturbatively measured values of P and R. For full QCD we use values 
extrap olated to t he chir al limit of the sea quark. Numerical values of P and R used in the calculation are given in 



Tables XXIV and XXV 



As an alternative we define the tadpole improved coupling constant with the usual procedure which only uses the 
plaquette P, 

IP 11 - -Nf 

0.2402 -I- 0.031497V/ -1- — — |— ^log(H. (C3) 



2 



The VWI quark mass is renormalized with 



niR^ Zrn[l+bm—] (C4) 
Mo/ Uo 



where 



gM^i) I 0.0400 - log(Ma)' ) , (C5) 
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and 

1 



6™ = --- 0.03235yM). (C6) 



For uq 

0.8412\^/'^ 



no = P^/^ = ( 1 - ^ j (C7) 



is used. 

The local pseudoscalar density P„ = 4'nl5''Pn is renormalized with 



f,f - 2kuoZp ( 1 + 6p— ) P„, (C8) 

"0/ 



where 



Zp = l + gl^ifi) (^-0.0523 + ^ log(Aia)2 j , (C9) 

and 

&P = l + 0.0397.g|g(A*). (CIO) 
The renormalized axial vector current Ajj, improved to 0{g^a), is obtained through 

A^^ = 2kuoZa (^1 + bA^^ (a„^ + c^a^Pn) , (Cll) 

where = ipj^ij^^/^tpn is the bare local current and 9^ the symmetric lattice derivative. Perturbative expressions 
for the renormalization factor and the improvement coefficients are 

Za = 1- 0.0215g|g(/i), (C12) 
&A = l + 0.0378g|g(/i), (C13) 
CA = -0.0038g|^(/.). (C14) 

Similarly, the renormalized vector current is obtained from the bare local vector current = V'n7/i'0n 
Tnf_Li, = i'J'^i-ii'i^n through 



= 2kuoZv ( 1 + — ) ( Km + cy 9.T„^. j . (C15) 



Here the perturbative results are 



Zv = l- 0m77g^{p), (C16) 



2 



bv = l + 0m82gj^{p), (C17) 



cv = -0.0097<?i^(M). (C18) 
Numerical values for coupling constants, Z factors and improvement coefficients are listed in Table pCXV^ 
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TABLE I. Overview of simulations. The scale a is fixed by Mp — 768.4 MeV from fit to vector mesons with Eq. 



p 


L '^ xT 


csw 


a [fm] 


La [fm] 






mps/my for sea 


quarks : 






1.80 
1.95 
2.10 
2.20 


12^x24 
16^x32 
24^x48 
24^x48 


1.60 
1.53 
1.47 
1.44 


0.2150(22) 
0.1555(17) 
0.1076(13) 
0.0865(33) 


2.580(26) 
2.489(27) 
2.583(31) 
2.076(79) 


0.807(1) 
0.804(1) 
0.806(1) 
0.799(3) 


6250 
7000 
4000 
2000 


0.753(1): 5000 
0.752(1): 7000 
0.755(2): 4000 
0.753(4): 2000 


0.694(2) 
0.690(1) 
0.691(3) 
0.705(5) 


:7000 
:7000 
:4000 
:2000 


0.547(4): 5250 
0.582(3): 5000 
0.576(3): 4000 
0.632(7): 2000 



TABLE IL Run parameters. The employed molecular dynamics (MD) integration schemes are introduced in Sec. II C, Ninv 
is the sum of iterations for inversions of and D in the evaluation of the fermionic force during HMC. A'Meas is the number 
of hadron propagator measurements. In brackets of A^Mcas the numbers of removed configurations are also given. 



/3 


Size 


csw 


K 


MD 


At 


Accept. 


Stop 


Ainv 


Hour/Traj. 


A^Traj 


A^Mcas 


1.80 


12^x24 


1.60 


0.1409 


c) 


0.033 


0.781 


10-1° 


64.8 


0.10 


6250 


1238(12) 




on 64 PU 




0.1430 


c) 


0.025 


0.807 


10-" 


87.2 


0.15 


5000 


990(10) 








0.1445 


c) 


0.0167 


0.840 


10-" 


119.5 


0.26 


3500 


690(10) 










a) 


0.0065 


0.809 


10-" 


120.4 


0.25 


3500 


692(8) 








0.1464 


a) 


0.0033 


0.764 


10-" 


263.6 


0.92 


4280 


839(17) 










b) 


0.0066 


0.714 


10-" 


256.9 


0.90 


970 


194(0) 


1.95 


16^x32 


1.53 


0.1375 


c) 


0.03125 


0.732 


10-" 


95.1 


0.10 


7000 


1400(0) 




on 256 PU 




0.1390 


c) 


0.025 


0.755 


10-" 


133.3 


0.15 


7000 


1395(5) 








0.1400 


c) 


0.0185 


0.761 


10-" 


187.4 


0.25 


7000 


1397(3) 








0.1410 


c) 


0.008 


0.820 


10-" 


331.8 


0.83 


5000 


1000(0) 


2.10 


24^x48 


1.47 


0.1357 


b) 


0.02 


0.759 


10-12 


151.3 


0.35 


4000 


798(2) 




on 512 PU 




0.1367 


b) 


0.016 


0.792 


10-12 


208.7 


0.57 


4000 


800(0) 








0.1374 


b) 


0.0143 


0.788 


10-12 


289.3 


0.82 


4000 


798(2) 








0.1382 


b) 


0.0075 


0.781 


10-12 


544.7 


2.72 


4000 


800(0) 


2.20 


24^x48 


1.44 


0.1351 


b) 


0.02 


0.758 


10-12 


192.0 


0.42 


2000 


400(0) 




on 512 PU 




0.1358 


b) 


0.016 


0.826 


10-12 


254.9 


0.67 


2000 


400(0) 








0.1363 


b) 


0.0143 


0.837 


10-12 


336.8 


0.94 


2000 


400(0) 








0.1368 


b) 


0.01 


0.859 


10-12 


505.6 


1.90 


2000 


400(0) 



TABLE III. Parameters A and B used for the smearing of quark sources. 

(3 = 1.80 P = 1.95 13 = 2.10 f3 = 2.20 

A TIB L25 r05 1.02 

B 0.50 0.50 0.17 0.125 



TABLE IV. Fit ranges for extraction of potential data, Eq. (|28|), and ranges of -Rmin and i?max used in potential fit, Eq. ( [29[ 
At /3 = 2.2, we use the same fit ranges as those taken ai l3 — 2.1. 





13 = 1.80 


P = 1.95 


(3 = 2.10 


t 


[2,4] 


[3,6] 


[4,8] 




[y2,x/2] 


[72, V6] 


[V5,3] 




[2^3,4] 


[3x/5,8] 


[9,6%/5] 
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TABLE V. String tension o and Sommer scale ro at simulated sea quark masses and in the chiral limit of the sea quark. k\ 
is the hopping parameter corresponding to the heaviest sea quark, to the lightest, and in the chiral limit mps = 
are obtained from extrapolations using Eqs. (H) and (|^). 



K, 


13^ 


1.80 


/3 = 


1.95 


fi = 


2.10 


fi = 


2.20 




a 




a 




a 




a 




Kl 


0.4115(96) 


1.716(35) 


0.2078(22) 


2.497(54) 


0.08949(99) 


3.843(16) 


0.05485(17) 


4.913(21) 


K2 


0.389(12) 


1.799(13) 


0.1859(29) 


2.651(42) 


0.07823(90) 


4.072(15) 


0.05107(26) 


5.073(19) 


«3 


0.3595(68) 


1.897(30) 


0.1633(23) 


2.821(29) 


0.07195(73) 


4.236(14) 


0.04760(31) 


5.237(22) 


K4 


0.3067(60) 


2.064(38) 


0.1436(25) 


3.014(33) 


0.06340(51) 


4.485(12) 


0.04474(23) 


5.410(21) 


'TIPS = 


0.2858(72) 


2.175(51) 


0.1295(25) 


3.210(52) 


0.05720(63) 


4.695(18) 


0.04072(29) 


5.656(33) 



TABLE VI. 


Estimate of intej 


grated autocorrelation times for the j 


^auge action S^, for the number 


of iterations for the 


inversion of the Dirac matrix Mnv, and for the effective pion mass m^.t 


fj at the onset of a plateau. 






K 


rs, [xl KMC Traj.] riv, 


„^ [xl HMC Traj.] r™ 


[x5 HMC Traj.] 


1.80 


0.1409 


4.2(4) 


3.4(3) 


0.7(1) 




0.1430 


5.5(1.0) 


4.6(7) 


0.7(1) 




0.1445 


6.9(9) 


5.7(8) 


0.7(2) 




0.1464 


9.2(1.6) 


7.8(9) 


0.9(2) 


1.95 


0.1375 


5.6(6) 


5.8(7) 


0.9(2) 




0.1390 


7.6(9) 


8.2(1.6) 


1.4(2) 




0.1400 


9.0(2.0) 


10.1(2.0) 


1.4(2) 




0.1410 


7.8(1.4) 


9.3(2.0) 


1.0(2) 


2.10 


0.1357 


3.8(5) 


6.7(1.5) 


1.5(3) 




0.1367 


4.0(7) 


9.4(2.4) 


1.5(3) 




0.1374 


3.1(4) 


8.1(2.1) 


1.0(1) 




0.1382 


5.4(1.0) 


11.0(2.1) 


1.7(4) 


2.20 


0.1351 


2.7(5) 


5.0(1.5) 


1.5(2) 




0.1358 


1.8(3) 


4.4(1.2) 


3.9(1.6) 




0.1363 


2.1(3) 


3.2(8) 


2.5(1.0) 




0.1368 


1.9(3) 


4.3(1.2) 


2.5(7) 



TABLE VII. Parameters of chiral fits to pseudoscalar meson masses as a function of Xjn with Eq. ( |40[ ) (first four rows) or 
as a function of the AWI quark mass with Eq. (t44) (last four rows). 



/3 






hs 






Cv 




C-vv 


1.80 


116/29 


0.147635(16) 


4.562(72) 


5.400(19) 


-11.51(38) 


-3.064(81) 


2.45(10) 


1.646(42) 


1.95 


26.6/29 


0.142065(13) 


2.655(69) 


4.169(12) 


-4.64(48) 


-0.846(72) 


4.379(99) 


1.333(42) 


2.10 


17.4/29 


0.138984(13) 


0.924(55) 


3.206(13) 


-1.40(49) 


0.96(12) 


4.38(17) 


1.121(91) 


2.20 


15.0/29 


0.137675(52) 


0.55(19) 


2.685(36) 


-1.9(1.9) 


1.79(38) 


4.52(44) 


1.04(27) 


a 












c' 


c' 




1.80 


75/33 






5.777(25) 




-1.335(90) 


1.99(15) 




1.95 


75/33 






4.393(28) 




-0.33(11) 


3.81(24) 




2.10 


57/33 






3.188(25) 




1.02(14) 


3.63(35) 




2.20 


32/33 






2.641(61) 




1.37(37) 


3.37(86) 
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TABLE VIII. Parameters of chiral fits to vector meson 


masses with Eq. (| 


48| (first four rows 


) or Eq. (^ (last four rows). 


a 






bX 


-fc>„ 






^ sv 


f.80 


25.5/30 


0.8241(85) 


0.206(15) 


4066f96l 


— 0517f87") 


-0.0193(48) 


-0.0471(66) 


1.95 


32.8/30 


0.5963(66) 


0.258(23) 


0.567(11) 


-0.072(21) 


-0.0443(91) 


-0.112(12) 


2.10 


43.0/30 


0.4124(51) 


0.327(34) 


0.907(24) 


-0.143(70) 


-0.265(40) 


-0.368(46) 


o on 


b.U/oU 


\j.66Z\l6 ) 


U.4D i 


1 0801^571 


— 40(40) 




— U. (K)\16 ) 


/3 






bX 








^ sv 


1.80 


22.0/30 


0.802(10) 


0.319(28) 


0.480(18) 


-0.120(20) 


-0.078(11) 


-0.086(12) 


1.95 


25.3/30 


0.5812(80) 


0.376(42) 


0.663(19) 


-0.131(37) 


-0.135(16) 


-0.155(17) 


2.10 


37.4/30 


0.4003(63) 


0.478(64) 


1.112(44) 


-0.192(88) 


-0.459(52) 


-0.350(45) 


2.20 


4.7/30 


0.320(16) 


0.70(26) 


1.32(11) 


-0.41(42) 


-0.59(15) 


-0.60(18) 


TABLE IX. Parameters of chiral fits to decuplet baryon masses with Eq. (|5l|). 


/3 


x'/A'df 




Bf 


Bv 


cf 




^ SV 


1.80 


13.5/46 


1.360(24) 


0.461(55) 


0.647(36) 


-0.116(33) 


-0.036(22) 


-0.090(25) 


1.95 


2.12/46 


1.036(17) 


0.384(65) 


0.816(26) 


-0.038(67) 


-0.034(25) 


-0.193(38) 


2.10 


7.82/46 


0.704(17) 


0.67(12) 


1.202(67) 


-0.52(23) 


-0.11(11) 


-0.48(13) 


2.20 


15.9/46 


0.527(28) 


1.20(34) 


1.64(14) 


-1.9(1.0) 


-0.75(34) 


-0.89(43) 






TABLE X. Parameters of chiral fits to octet bary( 


3n masses with Eq. (^). 




/3 


xVA'df 


A° 




^ V 




c? 


r° 






r<o 


C„ 




r° 

^ sv 


^ sv 


^ sv 


1.80 


28/72 


1.080(18) 


0.303(39) 


0.2945(96) 


-0.0685(96) 


-0.056(23) 


-0.0437(67) 






-0.0507(69) 


0.0027(41) 


0.0401(59) 


-0.0411(82) 


0.0225(70) 


0.0118(65) 


1.95 


17/72 


0.804(11) 


0.219(44) 


0.3799(69) 


-0.0959(62) 


0.009(45) 


-0.0855(77) 






-0.1163(78) 


0.0090(54) 


0.0952(57) 


-0.055(13) 


0.0334(73) 


0.0039(76) 


2.10 


59/72 


0.5418(83) 


0.376(64) 


0.576(13) 


-0.131(13) 


-0.17(13) 


-0.251(29) 






-0.389(29) 


0.051(24) 


0.289(22) 


-0.189(43) 


0.107(29) 


0.006(27) 


2.20 


12/72 


0.432(19) 


0.65(24) 


0.674(42) 


-0.170(36) 


-1.06(73) 


-0.46(14) 






-0.61(11) 


0.041(92) 


0.519(80) 


-0.32(17) 


0.20(14) 


0.03(12) 






TABLE XL 


Lattice spacings and hopping parameters and k.. 






P 




a [fm] 






Ka {Mk input) 




Ks {Mif, input) 


1.80 




0.2150(22) 


0.147540(16) 


0.143147(91) 




0.14192(16) 


1.95 




0.1555(17) 


0.141998(12) 


0.139279(59) 




0.138633(79) 


2.10 




0.1076(13) 


0.138933(12) 


0.137324(41) 




0.137105(61) 


2.20 




0.0865(33) 


0.137634(50) 


0.13642(11) 




0.13622(11) 
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TABLE XII. Meson and baryon masses at finite lattice spacings and in the continuum limit. Values in the continuum limit 
are obtained by a fit linear in the lattice spacing to data at /? = 1.8, 1.95 and 2.1. All masses are in GeV units. 



channel 


experiment 


13 = 1.8 


13 = 1.95 


P = 2.1 


P = 2.2 


Continuum 


iV 


0.9396 


1.016(16) 


1.040(13) 


1.016(17) 


1.007(43) 


1.034(36) 


A 


1.232 


1.270(23) 


1.332(18) 


1.310(30) 


1.225(56) 


1.392(58) 


Mk input 


Vss 




0.69154(21) 


0.69578(19) 


0.69769(29) 


0.69838(91) 




K* 


0.8961 


0.8685(16) 


0.8708(13) 


0.8813(19) 


0.8774(46) 


0.8902(38) 


<!> 


1.0194 


0.9660(27) 


0.9710(23) 


0.9895(33) 


0.9832(84) 


1.0066(67) 


A 


1.1157 


1.149(15) 


1.165(12) 


1.147(16) 


1.132(38) 


1.158(33) 


E 


1.1926 


1.183(14) 


1.202(11) 


1.183(15) 


1.169(38) 


1.197(32) 




1.3149 


1.295(13) 


1.304(11) 


1.292(14) 


1.271(35) 


1.298(30) 


E* 


1.3837 


1.376(20) 


1.431(16) 


1.411(27) 


1.336(53) 


1.485(51) 




1.5318 


1.481(18) 


1.529(15) 


1.512(25) 


1.443(51) 


1.577(47) 




1.6725 


1.583(17) 


1.627(15) 


1.612(23) 


1.548(50) 


1.671(44) 


input 


K 


0.4977 


0.5583(35) 


0.5506(28) 


0.5287(36) 


0.5355(98) 


0.5042(78) 


Vss 




0.7791(50) 


0.7738(41) 


0.7438(55) 


0.755(14) 




K* 


0.8961 


0.89607(50) 


0.89573(34) 


0.89698(35) 


0.89616(79) 


0.89778(86) 


A 


1.1157 


1.184(14) 


1.195(12) 


1.165(15) 


1.153(39) 


1.160(32) 


S 


1.1926 


1.225(13) 


1.239(11) 


1.205(15) 


1.195(39) 


1.202(30) 




1.3149 


1.367(13) 


1.365(11) 


1.329(13) 


1.314(37) 


1.302(28) 


E* 


1.3837 


1.406(19) 


1.455(16) 


1.426(26) 


1.355(52) 


1.488(49) 




1.5318 


1.538(17) 


1.577(15) 


1.541(23) 


1.480(50) 


1.583(44) 




1.6725 


1.666(16) 


1.698(14) 


1.654(21) 


1.601(49) 


1.680(41) 



TABLE XIII. Parameters of quenched QCD simulations. Coupling constants /3 are chosen, so that measured values of a 
correspond to the ones in full QCD given in Table vA 



p 


Size 


csw 


a [fm] 


La [fm] 




r(3 


2.187 


16^x32 


1.439 


0.2004(20) 


3.206(31) 


0.2157(32) 


2.494(35) 


2.214 


16^x32 


1.431 


0.1903(19) 


3.045(31) 


0.1949(25) 


2.621(43) 


2.247 


16^x32 


1.422 


0.1807(18) 


2.891(29) 


0.1713(18) 


2.801(28) 


2.281 


16^x32 


1.412 


0.1765(20) 


2.824(32) 


0.1487(17) 


3.001(36) 


2.334 


16^x32 


1.398 


0.1632(16) 


2.611(26) 


0.1241(14) 


3.289(23) 


2.416 


24^x48 


1.378 


0.1446(18) 


3.471(42) 


0.0921(10) 


3.824(13) 


2.456 


24^x48 


1.370 


0.1328(13) 


3.188(30) 


0.0800(16) 


4.080(16) 


2.487 


24^x48 


1.363 


0.1284(14) 


3.081(34) 


0.0725(11) 


4.286(15) 


2.528 


24^x48 


1.355 


0.1206(13) 


2.895(30) 


0.0637(11) 


4.570(21) 


2.575 


24^x48 


1.345 


0.1130(11) 


2.713(27) 


0.0561(7) 


4.887(16) 



TABLE XIV. Parameters of chiral fits to meson masses in quenched QCD with Eqs. (BJ, (BS) and (M). 



p 




bv 


Cv 




b'v 


c' 






2.187 


0.141666(12) 


4.660(25) 


-1.59(13) 


1.631(97) 


4.684(48) 


0.66(32) 


0.7735(77) 


0.4229(87) 


2.214 


0.140999(15) 


4.496(28) 


-1.42(13) 


1.612(84) 


4.582(44) 


0.20(26) 


0.7349(74) 


0.447(10) 


2.247 


0.140239(19) 


4.408(40) 


-1.68(24) 


1.63(14) 


4.370(44) 


0.25(28) 


0.6975(71) 


0.480(11) 


2.281 


0.139587(15) 


4.211(37) 


-1.24(25) 


1.50(14) 


4.162(36) 


0.38(27) 


0.6816(77) 


0.473(12) 


2.334 


0.138728(13) 


3.849(27) 


-0.44(18) 


1.37(12) 


3.854(42) 


0.72(29) 


0.6302(63) 


0.513(11) 


2.416 


0.137633(7) 


3.434(17) 


0.32(20) 


1.21(14) 


3.388(34) 


1.22(30) 


0.5586(69) 


0.571(16) 


2.456 


0.137179(6) 


3.258(15) 


0.70(14) 


1.05(11) 


3.220(26) 


1.06(29) 


0.5128(49) 


0.647(13) 


2.487 


0.136852(7) 


3.168(20) 


0.51(16) 


1.07(11) 


3.094(22) 


1.20(24) 


0.4956(56) 


0.662(16) 


2.528 


0.136493(7) 


2.951(23) 


1.33(16) 


0.81(12) 


2.918(20) 


1.05(21) 


0.4656(49) 


0.708(15) 


2.575 


0.136116(8) 


2.781(18) 


1.60(18) 


0.69(14) 


2.776(28) 


0.81(30) 


0.4364(44) 


0.757(16) 
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TABLE XV. Meson and baryon masses at finite lattice spacings and in the continuum limit in quenched QCD. All masses 
are in GeV units. 

Mk input 



niK' 



niN 



0.8690(11 
0.8694(14 
0.8710(15 
0.8675(16 
0.8684(1.3 
0.8673(15 
0.8712(12 
0.8699(14 
0.8706(12 
0.8709(12 
0.8728(21 



0.9695(23) 
0.9704(28) 
0.9736(31) 
0.9665(32) 
0.9684(25) 
0.9662(31) 
0.9740(25) 
0.9714(28) 
0.9729(23) 
0.9733(24) 
0.9773(42) 



0.991(18) 
0.987(17) 
1.014(15) 
0.958(19) 
0.953(17) 
0.943(17) 
0.961(15) 
0.925(15) 
0.983(18) 
0.943(15) 
0.873(28) 



1.173(14) 
1.169(13) 
1.189(11) 
1.152(13) 
1.143(13) 
1.133(12) 
1.152(11) 
1.119(12) 
1.156(14) 
1.140(12) 
1.079(21) 



1.126(14) 
1.133(12) 
1.135(12) 
1.092(14) 
1.096(12) 
1.076(13) 
1.097(12) 
1.069(11) 
1.112(13) 
1.091(11) 
1.024(20) 



1.277(10) 

1.2893(90) 

1.2802(89) 

1.248(11) 

1.2561(88) 

1.2322(99) 

1.2547(89) 

1.2305(87) 

1.2587(94) 

1.2545(80) 

1.196(15) 



1.324(18) 
1.339(19) 
1.334(16) 
1.274(19) 
1.305(15) 
1.283(20) 
1.302(20) 
1.253(20) 
1.298(17) 
1.289(17) 
1.219(30) 



1.420(15) 
1.434(16) 
1.430(14) 
1.379(16) 
1.403(14) 
1.384(17) 
1.402(16) 
1.360(16) 
1.398(14) 
1.391(14) 
1.331(25) 



1.516(12) 
1.529(14) 
1.526(11) 
1.483(13) 
1.501(12) 
1.486(14) 
1.503(13) 
1.467(13) 
1.499(12) 
1.493(11) 
1.443(21) 



1.612(10) 
1.625(11) 
1.623(10) 
1.587(11) 
1.599(10) 
1.587(12) 
1.603(11) 
1.574(11) 
1.599(10) 
1.595(10) 
1.555(17) 



input 



rriK 



ruN 



rriT, 



rriA 



mn 



0.5507(29 
0.5496(35 
0.5458(37 
0.5545(41 
0.5518(32 
0.5544(39 
0.5448(30 
0.5481(34; 
0.5462(28 
0.5456(29 
0.5400(52 



0.89100(25) 
0.89090(23) 
0.89145(32) 
0.89097(34) 
0.89013(27) 
0.88957(27) 
0.88991(22) 
0.89011(30) 
0.88974(33) 
0.88978(31) 
0.88760(48) 



0.991(18) 
0.987(17) 
1.014(15) 
0.958(19) 
0.953(17) 
0.943(17) 
0.961(15) 
0.925(15) 
0.983(18) 
0.943(15) 
0.873(28) 



1.212(13) 
1.207(12) 
1.223(11) 
1.195(11) 
1.184(12) 
1.176(10) 
1.188(10) 
1.158(10) 
1.190(13) 
1.177(11) 
1.113(20) 



1.158(13) 

1.165(11) 

1.163(11) 

1.128(13) 

1.130(11) 

1.111(11) 

1.126(11) 

1.101(9) 

1.140(12) 

1.123(10) 

1.052(18) 



1.340(8) 
1.353(8) 
1.336(8) 
1.319(9) 
1.323(7) 
1.303(6) 
1.313(7) 
1.294(6) 
1.317(8) 
1.316(7) 
1.250(12) 



1.324(18) 
1.339(19) 
1.334(16) 
1.274(19) 
1.305(15) 
1.283(20) 
1.302(20) 
1.253(20) 
1.298(17) 
1.289(17) 
1.219(30) 



1.444(13) 
1.457(15) 
1.452(13) 
1.406(14) 
1.428(12) 
1.412(15) 
1.424(15) 
1.385(15) 
1.421(13) 
1.414(13) 
1.355(23) 



1.564(10) 

1.576(11) 

1.569(10) 

1.539(10) 

1.551(9) 

1.540(10) 

1.547(11) 

1.517(10) 

1.544(10) 

1.539(9) 

1.490(17) 



1.683(8) 
1.694(8) 
1.687(9) 
1.671(8) 
1.673(7) 
1.669(6) 
1.669(9) 
1.650(7) 
1.668(7) 
1.664(8) 
1.622(13) 



TABLE XVL Renormalized quark masses (in MeV) in the MS-scheme at /i = 2 GcV at finite lattice spacings in full and 
quenched QCD. Values in the continuum limit obtained with separate linear fits to cacli definition arc also listed. For full QCD 
data at i = 2.2 wt^rt^ not iuciudt'd in these fits. 

7) _VWI _VWI,PQ IIaWI IIvwrTT^A ™VWI,PQ / _AWI / _VWI i i\ _VWI,PQ f i\ _AWI i i\ 



Nf = 2 Full QCD 

1.80 2.277(27) 4.183(42) 3.322(37) 102.92(92) 104.54(93) 88.0(1.0) 129.1(2.2) 130.7(2.2) 113.9(2.4) 

1.95 2.489(38) 4.064(43) 3.321(38) 100.65(98) 102.08(99) 87.2(1.0) 123.1(1.7) 124.5(1.7) 109.8(1.7) 

2.10 2.966(55) 3.816(47) 3.344(46) 95.6(1.1) 96.4(1.1) 87.0(1.2) 108.0(2.2) 108.8(2.2) 100.0(2.2) 

2.20 3.11(22) 3.75(15) 3.35(15) 94.4(3.5) 95.0(3.5) 86.9(3.9) 109.4(4.7) 110.0(4.7) 102.6(5.0) 

o-> 3.47(10) 3.50(10) 3.36(9) 89.4(2.3) 89.5(2.3) 85.8(2.4) 90.1(4.9) 90.3(4.9) 88.1(4.9) 

X^/NuF 10.8 2^4 007 2^1 2^7 003 6^0 05 2.4 

Nf = Quenched QCD 

2.187 4.429(50) - 3.873(53) 109.8(1.2) - 100.7(1.2) 133.5(2.5) - 125.7(2.6) 

2.214 4.387(47) - 3.791(52) 109.1(1.1) - 99.1(1.2) 132.2(2.6) - 124.1(2.9) 

2.247 4.273(59) - 3.802(53) 107.0(1.3) - 99.3(1.2) 128.2(2.7) - 122.0(2.8) 

2.281 4.374(63) - 3.913(52) 109.2(1.4) - 102.0(1.2) 134.8(3.2) - 129.8(3.3) 

2.334 4.458(47) - 3.950(56) 110.9(1.0) - 102.6(1.3) 135.5(2.5) - 129.4(2.7) 

2.416 4.481(57) - 4.045(60) 111.4(1.3) - 104.5(1.3) 137.3(3.3) - 132.6(3.4) 

2.456 4.378(45) - 3.955(43) 109.9(1.0) - 102.4(1.0) 1302(2.3) - 125.5(2.4) 

2.487 4.363(56) - 3.994(51) 109.1(1.3) - 103.3(1.2) 131.8(2.8) - 127.7(2.9) 

2.528 4.426(53) - 4.013(47) 110.1(1.2) - 103.9(1.1) 132.0(2.3) - 128.2(2.5) 

2.575 4.425(53) - 3.984(53) 110.2(1.1) - 103.5(1.2) 131.8(2.4) - 127.8(2.5) 

a-i-O 4.449(87) - 4.269(86) 111.2(1.9) - 109.4(2.0) 130.8(4.2) - 132.4(4.4) 

X^/NoF 1.2 - 1.0 1.0 - 0.9 1.0 - 1.0 
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TABLE XVII. Breakdown of contributions to total error of full QCD quark masses in the continuum limit. 



Stat. 



Cont. Ext. 



Chiral 



Z factor 



rus (Mk input) 
rris {M^ input) 



+2.6% 
-2.6% 
+2.4% 
-2.4% 
+4.8% 
-4.8% 



+ 1.7% 
-2.3% 
+ 1.4% 
-2.8% 
+0.9% 
-1.6% 



+ 1.2% 
-2.3% 
+1.6% 
-2.2% 
+ 1.5% 
-7.6% 



+2.3% 
-5.0% 
+2.2% 
-5.6% 
+1.7% 
-6.9% 



TABLE XVIIL Continuum limit quark masses in the MS scheme at /i = 2 GeV (in MeV) . 







action 


rriud 


Mk input 


rus 


input 




= 2 


impr. 


o 44+0.14 


88^^ 






Nf 


= 


impr. 


^"-""-O.lT 


iwtl 




132t* 


Nf 


= 


stand. iD 


4.57(18) 


116(3) 




144(6) 



TABLE XIX. Parameters of chiral fits to pseudoscalar decay constants (upper part) and vector decay constants (lower part) 
with Eq. 



p 


XV^DF 




Bf 




c! 


CiT 


^ SV 


1.80 


47.1/30 


0.2082(52) 


0.086(13) 


0.1309(61) 


-0.0228(84) 


-0.0110(32) 


-0.0069(41) 


1.95 


14.6/30 


0.1207(54) 


0.087(22) 


0.1502(56) 


-0.029(21) 


-0.0253(44) 


-0.0266(68) 


2.10 


14.8/30 


0.0696(37) 


0.168(35) 


0.189(12) 


-0.218(73) 


-0.101(23) 


-0.029(30) 


2.20 


10.6/30 


0.051(12) 


0.29(15) 


0.211(32) 


-0.79(44) 


-0.123(92) 


-0.08(15) 


1.80 


11.0/30 


0.362(10) 


0.226(24) 


0.065(13) 


-0.072(15) 


0.009(7) 


-0.016(10) 


1.95 


17.3/30 


0.2105(73) 


0.184(28) 


0.061(12) 


-0.068(27) 


0.020(11) 


-0.015(15) 


2.10 


8.0/30 


0.1290(51) 


0.160(41) 


0.130(22) 


-0.092(90) 


-0.049(39) 


-0.080(52) 


2.20 


9.7/30 


0.0970(90) 


0.25(10) 


0.128(45) 


-0.29(33) 


0.02(12) 


-0.24(16) 



TABLE XX. Decay constants at finite lattice spacings in full and quenched QCD. All decay constants are in GeV units. 





.U 


fK (K) 


fK (0) 


Fp 


Fk' (K) 


Fk' (0) 


F^ (K) 


F^ (0) 


Nf = 2 Full QCD 


1.80 


0.1954(51) 


0.2273(45) 


0.2359(43) 


0.3378(66) 


0.3546(61) 


0.3595(61) 


0.3726(65) 


0.3833(65) 


1.95 


0.1565(70) 


0.1832(65) 


0.1896(64) 


0.2705(73) 


0.2823(64) 


0.2854(63) 


0.2952(60) 


0.3019(58) 


2.10 


0.1311(66) 


0.1542(63) 


0.1573(63) 


0.2394(73) 


0.2555(65) 


0.2577(64) 


0.2708(65) 


0.2749(65) 


2.20 


0.120(26) 


0.141(26) 


0.145(26) 


0.224(16) 


0.237(14) 


0.2397(14) 


0.251(14) 


0.2555(14) 








Nf 


= Quenched QCD 








2.187 


0.1695(44) 


0.1912(36) 


0.1966(33) 


0.2861(44) 


0.3029(32) 


0.3070(30) 


0.3197(26) 


0.3280(25) 


2.214 


0.1622(39) 


0.1841(32) 


0.1894(30) 


0.2761(38) 


0.2917(29) 


0.2955(28) 


0.3074(30) 


0.3149(31) 


2.247 


0.1574(42) 


0.1797(35) 


0.1847(27) 


0.2706(37) 


0.2876(30) 


0.2914(29) 


0.3046(33) 


0.3122(32) 


2.281 


0.1477(34) 


0.1722(29) 


0.1787(27) 


0.2704(38) 


0.2834(29) 


0.2868(27) 


0.2963(27) 


0.3033(26) 


2.334 


0.1511(43) 


0.1716(37) 


0.1768(34) 


0.2601(30) 


0.2713(22) 


0.2742(22) 


0.2825(22) 


0.2882(23) 


2.416 


0.1407(40) 


0.1607(31) 


0.1661(28) 


0.2471(54) 


0.2557(37) 


0.2581(34) 


0.2644(28) 


0.2690(29) 


2.456 


0.1482(40) 


0.1661(34) 


0.1700(33) 


0.2332(44) 


0.2460(34) 


0.2488(32) 


0.2588(31) 


0.2645(31) 


2.487 


0.1391(37) 


0.1586(31) 


0.1632(29) 


0.2467(42) 


0.2558(30) 


0.2579(29) 


0.2648(27) 


0.2691(28) 


2.528 


0.1436(48) 


0.1626(40) 


0.1669(38) 


0.2293(45) 


0.2422(34) 


0.2451(32) 


0.2551(27) 


0.2610(25) 


2.575 


0.1476(55) 


0.1658(43) 


0.1699(40) 


0.2417(37) 


0.2487(29) 


0.2503(27) 


0.2557(26) 


0.2589(27) 
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TABLE XXI. Meson masses and AWI quark masses. 



'*val 


(2) 
'*val 


mps 


[imin, imax] 


X^/NdF mv [tminjtmax] / Nu¥ 


mps /mv 




13 = 1.80, Ksea = 0.1409 



0.1409 


0.1409 


1.15601(61) 


[5,12] 


1.1(9) 


1.4330(13) 


[5,11] 


1.4 


1.2) 


0.80669(71) 


0.22483 


'43) 


0.1409 


0.1430 


1.09336(63) 


[5,12] 


1.2(9) 


1..3908(15) 


[5,11] 


1.4 


1.1) 


0.78616(80) 


0.20056 


'40) 


0.1409 


0.1445 


1.04660(64) 


[5,12] 


1.3(9) 


1.3603(16) 


[5,11] 


1.3 


1.1) 


0.76938(88) 


0.18328 


'39) 


0.1409 


0.1464 


0.98441(68) 


[5,12] 


1.5(1.0) 


1.3217(18) 


[5,11] 


1.0 


1.0) 


0.7448(10) 


0.16141 


;38) 


0.1409 


0.1474 


0.94996(71) 


[5,12] 


1.6(1.1) 


1.3016(20) 


[5,11] 


0.7 


8) 


0.7298(11) 


0.14984 


;38) 


0.1430 


0.1430 


1.02733(65) 


[5,12] 


1.3(1.0) 


1.3479(16) 


[5,11] 


1.3 


1.1) 


0.76219(91) 


0.17649 


;38) 


0.1445 


0.1445 


0.92555(69) 


[5,12] 


1.5(1.1) 


1.2856(20) 


[5,11] 


1.3 


1.0) 


0.7199(11) 


0.14231 


'35) 


0.1464 


0.1464 


0.77767(78) 


[5,12] 


1.7(1.1) 


1.2051(29) 


[5,11] 


1.0 


9) 


0.6453(16) 


0.09916 


;37) 


0.1474 


0.1474 


0.6843(11) 


[5,12] 


0.9(8) 


1.1627(45) 


[5,11] 


0.7 


8) 


0.5885(24) 


0.07564 


;5o) 


3 = 1.80, Ksea = 0.1430 


0.1409 


0.1409 


1.11574(82) 


[6,12] 


1.3(1.1) 


1.3930(19) 


[6,12] 


1.1 


1.0) 


0.80091(97) 


0.21272 


;7i) 


0.1409 


0.1430 


1.05106(85) 


[6,12] 


1.2(1.0) 


1.3497(21) 


[6,12] 


1.0 


9) 


0.7787(11) 


0.18831 


;66) 


0.1430 


0.1430 


0.98267(89) 


[6,12] 


1.0(9) 


1.3057(24) 


[6,12] 


0.8 


0.8) 


0.7526(13) 


0.16412 


;6i) 


0.1430 


0.1445 


0.93112(92) 


[6,12] 


1.0(9) 


1.2743(27) 


[6,12] 


0.7 


8) 


0.7307(14) 


0.14692 


;57) 


0.1430 


0.1464 


0.8616(10) 


[6,12] 


1.1(9) 


1.2348(33) 


[6,12] 


0.6 


7) 


0.6978(18) 


0.12517 


;52) 


0.1430 


0.1474 


0.8225(11) 


[6,12] 


1.2(1.0) 


1.2149(40) 


[6,12] 


0.7 


8) 


0.6770(22) 


0.11359 


;5i) 


0.1445 


0.1445 


0.87676(96) 


[6,12] 


1.0(9) 


1.2424(31) 


[6,12] 


0.5 


7) 


0.7057(17) 


0.12983 


'53) 


0.1464 


0.1464 


0.7204(11) 


[5.12] 


1.4(1.1) 


1.1588(42) 


[5.11] 


0.4 


6) 


0.6217(23) 


0.08687 


54) 


0.1171 


0.1171 


0.()201(17) 


[•-',12] 


2.()(1.1) 


l.lir)()(()8) 


[5.11] 


0.1 


1) 


().5.55cS(:i7) 


().()():i()5 


78) 


(3 = 1.80, Ksoa = 0.1445 


0.1409 


0.1409 


1.07014(71) 


[6,12] 


2.1(1.3) 


1.3415(16) 


[6,12] 


0.4 


6) 


0.79774(85) 


0.19874 


;77) 


0.1409 


0.1445 


0.95358(75) 


[6,12] 


1.8(1.3) 


1.2637(20) 


[6,12] 


0.4 


6) 


0.7546(11) 


0.15664 


;63) 


0.1430 


0.1430 


0.93270(76) 


[6,12] 


1.8(1.3) 


1.2502(21) 


[6,12] 


0.4 


6) 


0.7460(12) 


0.14976 


;62) 


0.1430 


0.1445 


0.87921(78) 


[6,12] 


1.7(1.2) 


1.2174(23) 


[6,12] 


0.4 


6) 


0.7222(14) 


0.13250 


;56) 


0.1445 


0.1445 


0.82249(82) 


[6,12] 


1.6(1.2) 


1.1844(27) 


[6,12] 


0.5 


7) 


0.6945(16) 


0.11517 


'40) 


0.1445 


0.1464 


0.74507(83) 


[5,12] 


1.5(1.1) 


1.1433(35) 


[6,12] 


0.7 


8) 


0.6517(20) 


0.09359 


37) 


0.1445 


0.1474 


0.69993(92) 


[5,12] 


1.5(1.1) 


1.1227(43) 


[6,12] 


0.7 


8) 


0.6234(25) 


0.08213 


;38) 


0.1464 


0.1464 


0.65780(95) 


[5,12] 


1.7(1.2) 


1.1021(37) 


[5,12] 


0.5 


6) 


0.5969(21) 


0.07237 


;34) 


0.1474 


0.1474 


0.5464(16) 


[5,12] 


3.0(1.4) 


1.0600(65) 


[5,12] 


0.5 


7) 


0.5154(34) 


0.04904 


;46) 


13 = 1.80, Ksoa = 0.1464 


0.1409 


0.1409 


0.9818(11) 


[6,12] 


1.8(1.3) 


1.2427(20) 


[6,12] 


1.2 


1.0) 


0.7901(11) 


0.17153 


;77) 


0.1409 


0.1464 


0.7873(12) 


[6,12] 


1.1(1.0) 


1.1157(40) 


[6,12] 


1.5 


1.2) 


0.7056(26) 


0.10892 


;57) 


0.1430 


0.1430 


0.8346(12) 


[6,12] 


1.3(1.1) 


1.1436(27) 


[6,12] 


1.1 


1.0) 


0.7298(17) 


0.12322 


'64) 


0.1430 


0.1464 


0.6993(13) 


[6,12] 


0.9(9) 


1.0651(48) 


[6,12] 


1.7 


1.3) 


0.6565(31) 


0.08532 


'45) 


0.1445 


0.1445 


0.7152(13) 


[6,12] 


0.9(9) 


1.0725(38) 


[6,12] 


1.4 


1.1) 


0.6669(24) 


0.08951 


'46) 


0.1445 


0.1464 


0.6300(14) 


[6,12] 


0.8(8) 


1.0286(59) 


[6,12] 


1.7 


1.3) 


0.6125(37) 


0.06884 


;4o) 


0.1464 


0.1464 


0.5306(17) 


[6,12] 


0.7(8) 


0.9708(71) 


[5,11] 


1.5 


1.2) 


0.5466(44) 


0.04822 


;36) 


0.1464 


0.1474 


0.4666(25) 


[6,12] 


0.7(8) 


0.944(10) 


[5,12] 


1.3 


1.1) 


0.4940(64) 


0.03690 


'40) 


0.1474 


0.1474 


0.3872(52) 


[6,12] 


1.1(1.1) 


0.9307(81) 


[4,9] 


0.7 


8) 


0.4161(70) 


0.02478 


;63) 


(3 = 1.95, fvsea = 0.1375 


0.1375 


0.1375 


0.89400(52) 


[7,16] 


1.8(1.0) 


1.1113(13) 


[7,16] 


1.2 


8) 


0.80446(80) 


0.16112 


'63) 


0.1375 


0.1390 


0.83986(54) 


[7,16] 


1.8(1.0) 


1.0728(14) 


[7,16] 


0.9 


7) 


0.78289(92) 


0.14242 


58) 


0.1375 


0.1400 


0.80222(56) 


[7,16] 


1.9(1.0) 


1.0470(15) 


[7,16] 


0.7 


6) 


0.7662(10) 


0.12998 


;55) 


0.1375 


0.1410 


0.76315(58) 


[7,16] 


1.9(1.0) 


1.0213(17) 


[7,16] 


0.6 


6) 


0.7472(12) 


0.11755 


;53) 


0.1375 


0.1415 


0.74298(59) 


[7,16] 


1.8(1.0) 


1.0086(18) 


[7,16] 


0.7 


7) 


0.7366(13) 


0.11132 


;52) 


0.1390 


0.1390 


0.78290(56) 


[7,16] 


1.9(1.0) 


1.0337(16) 


[7,16] 


0.6 


6) 


0.7574(11) 


0.12387 


;53) 


0.1400 


0.1400 


0.70121(59) 


[7,16] 


1.8(1.0) 


0.9810(19) 


[7,16] 


0.5 


5) 


0.7148(14) 


0.09928 


'47) 


0.1410 


0.1410 


0.61020(63) 


[7,16] 


1.6(1.0) 


0.9274(20) 


[6,16] 


0.7 


6) 


0.6580(15) 


0.07484 


'42) 


0.1415 


0.1415 


0.55935(67) 


[7,16] 


1.3(9) 


0.9010(24) 


[6,16] 


0.8 


7) 


0.6208(18) 


0.06264 


;4i) 


(3 = 1.95, Ksea = 0.1390 


0.1375 


0.1375 


0.84401(65) 


[7,16] 


1.6(1.0) 


1.0495(15) 


[7,14] 


0.8 


8) 


0.80423(99) 


0.14824 


;57) 


0.1375 


0.1390 


0.78790(66) 


[7,16] 


1.6(1.0) 


1.0092(17) 


[7,14] 


0.7 


7) 


0.7807(11) 


0.12948 


'52) 


0.1390 


0.1390 


0.72857(68) 


[7,16] 


1.5(9) 


0.9686(18) 


[7.14] 


0.5 


6) 


0.7522(12) 


0.11090 


'46) 


0.1390 


0.1400 


0.68683(69) 


[7,16] 


1.3(9) 


0.9415(19) 


[7,14] 


0.5 


6) 


0.7295(13) 


0.09855 


;4i) 


0.1390 


0.1410 


0.64291(72) 


[7,16] 


1.2(8) 


0.9146(21) 


[7,14] 


0.5 


6) 


0.7029(15) 


0.08621 


;37) 


0.1390 


0.1415 


0.61988(74) 


[7,16] 


1.1(8) 


0.9015(23) 


[7,14] 


0.5 


6) 


0.6876(16) 


0.08002 


;35) 
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0.1400 


0.1400 


0.64284(71) 


[7,16] 


1.2(8) 


0.9141 


20) 


[7,14] 


0.4(6) 


0.7032(15) 


0.08628(37) 


0.1410 


0.1410 


0.54556(78) 


[7,16] 


1.0(8) 


0.8597 


27) 


[7,14] 


0.7(7) 


0.6346(20) 


0.06180(30) 


0.1415 


0.1415 


0.48957(89) 


[7,16] 


1.1(8) 


0.8325 


36) 


[7,14] 


1.1(9) 


0.5881(25) 


0.04954(28) 










13 = 


1.95, Ksea ~ 


= 0.1400 










0.1375 


0.1375 


0.80471(59) 


[6,16] 


1.2(7) 


1.0020 


13) 


[6,16] 


1.1(7) 


0.80308(86) 


0.13868(56) 


0.1375 


0.1400 


0.70643(63) 


[6,16] 


0.9(6) 


0.9327 


16) 


[6,16] 


1.1(7) 


0.7574(11) 


0.10705(46) 


0.1390 


0.1390 


0.68539(63) 


[6,16] 


0.8(6) 


0.9184 


16) 


[6,16] 


1.1(7) 


0.7463(12) 


0.10087(44) 


0.1390 


0.1400 


0.64192(66) 


[6,16] 


0.7(6) 


0.8906 


18) 


[6,16] 


1.1(8) 


0.7208(13) 


0.08839(40) 


0.1400 


0.1400 


0.59580(69) 


[6,16] 


0.7(6) 


0.8630 


20) 


[6,13] 


0.7(7) 


0.6904(14) 


0.07602(36) 


0.1400 


0.1410 


0.54639(74) 


[6,16] 


0.6(5) 


0.8354 


23) 


[6,13] 


1.2(8) 


0.6540(17) 


0.06369(33) 


0.1400 


0.1415 


0.51994(80) 


[6,16] 


0.6(5) 


0.8224 


25) 


[6,16] 


1.2(8) 


0.6322(19) 


0.05750(32) 


0.1410 


0.1410 


0.49232(82) 


[6,16] 


0.6(6) 


0.8082 


26) 


[6.16] 


1.1(9) 


0.6001(20) 


0.05143(30) 


0.111") 


0.1115 


0. l:ill(10) 


[(U(i] 


().8(()) 


0.7820 


3:!) 


[6^i:5] 


1.1(1.1) 


()..5.512(2()) 


().():!9()()(29) 










/? = 


1.95, /^sca — 


= 0.1410 










0.1375 


0.1375 


0.75717(73) 


[7,16] 


0.3(4) 


0.9416 


16) 


[6,15] 


1.3(9) 


0.8041(12) 


0.12759(44) 


0.1375 


0.1410 


0.61114(83) 


[6,16] 


0.3(4) 


0.8393 


24) 


[6,14] 


2.0(1.2) 


0.7281(20) 


0.08329(49) 


0.1390 


0.1390 


0.63303(82) 


[7,16] 


0.3(4) 


0.8527 


21) 


[6,14] 


1.9(1.2) 


0.7424(18) 


0.08981(39) 


0.1390 


0.1410 


0.53898(87) 


[6,16] 


0.4(5) 


0.7938 


30) 


[6,14] 


2.2(1.2) 


0.6790(25) 


0.06473(44) 


0.1400 


0.1400 


0.53870(85) 


[6,16] 


0.4(5) 


0.7929 


29) 


[6,14] 


2.3(1.2) 


0.6794(24) 


0.06477(44) 


0.1400 


0.1410 


0.48589(91) 


[6,16] 


0.6(6) 


0.7636 


34) 


[6,14] 


2.1(1.2) 


0.6363(28) 


0.05243(41) 


0.1410 


0.1410 


0.42700(98) 


[6,16] 


0.8(8) 


0.7339 


40) 


[6,14] 


1.7(1.1) 


0.5819(32) 


0.04020(37) 


0.1410 


0.1415 


0.3942(10) 


[6,16] 


1.1(8) 


0.7191 


43) 


[6,14] 


1.3(1.0) 


0.5481(34) 


0.03406(35) 


0.1415 


0.1415 


0.3582(11) 


[6,16] 


1.6(1. 


0) 0.7040 


48) 


[6,14] 


1.1(9) 


0.5089(37) 


0.02793(34) 










/3 = 


2.10, Ksca — 


- 0.1357 










0.1357 


0.1357 


0.63010(61) 


[10,24] 


0.9(7) 


0.7822 


16) 


[10,24] 


1.1(7) 


0.8055(14) 


0.10748(51) 


0.1357 


0.1367 


0.58676(62) 


[10,24] 


1.0(7) 


0.7509 


'12) 


[10,24] 


1.2(7) 


0.7814(16) 


0.09386(48) 


0.1357 


0.1374 


0.55502(64) 


[10,24] 


1.0(7) 


0.7292 


19) 


[10,24] 


1.3(7) 


0.7611(17) 


0.08407(44) 


0.1357 


0.1382 


0.51712(56) 


[10,24] 


1.2(7) 


0.7056 


21) 


[10,24] 


1.4(8) 


0.7329(20) 


0.07311(43) 


0.1357 


0.1385 


0.50240(70) 


[10,24] 


1.3(8) 


0.6974 


22) 


[10,24] 


1.5(8) 


0.7204(22) 


0.06898(43) 


0.1367 


0.1367 


0.54107(63) 


[10,24] 


1.0(7) 


0.7194 


19) 


[10,24] 


1.2(7) 


0.7521(18) 


0.08031(44) 


0.1374 


0.1374 


0.47157(64) 


[10,24] 


1.2(7) 


0.6762 


'23) 


[10,24] 


1.5(8) 


0.6974(22) 


0.06119(37) 


0.1382 


0.1382 


0.37964(71) 


[10,24] 


1.4(9) 


0.6273 


28) 


[9,21] 


1.7(1.0) 


0.6052(29) 


0.03957(29) 


0.1385 


0.1385 


0.33926(75) 


[10,24] 


1.5(9) 


0.6092 


38) 


[9,21] 


1.9(1.0) 


0.5569(38) 


0.03144(24) 










P = 


2.10, Ksea — 


= 0.1367 










0.1357 


0.1357 


0.60740(64) 


[10,24] 


1.6(9) 


0.7508 


13) 


[10,24] 


1.8(9) 


0.8090(13) 


0.10267(44) 


0.1357 


0.1367 


0.56332(66) 


[10,24] 


1.6(9) 


0.7179 


15) 


[10,24] 


1.6(9) 


0.7846(16) 


0.08912(41) 


0.1367 


0.1367 


0.51671(67) 


[10,24] 


1.5(9) 


0.6843 


16) 


[9,24] 


1.5(9) 


0.7551(17) 


0.07564(38) 


0.1367 


0.1374 


0.48208(68) 


[10,24] 


1.4(9) 


0.6612 


19) 


[9,24] 


1.6(8) 


0.7291(20) 


0.06617(35) 


0.1367 


0.1382 


0.44003(70) 


[10,24] 


1.2(8) 


0.6352 


24) 


[8,24] 


1.9(9) 


0.6928(25) 


0.05525(31) 


0.1367 


0.1385 


0.42339(71) 


[10,24] 


1.1(7) 


0.6262 


26) 


[8,24] 


1.9(9) 


0.6761(27) 


0.05114(29) 


0.1374 


0.1374 


0.44539(68) 


[10,24] 


1.3(8) 


0.6373 


23) 


[8,24] 


2.0(9) 


0.6989(24) 


0.05671(32) 


0.1382 


0.1382 


0.34991(68) 


[10,24] 


0.8(7) 


0.5801 


23) 


[7,16] 


2.4(1.4) 


0.6032(25) 


0.03476(24) 


0.1385 


0.1385 


0.30689(70) 


[10,24] 


0.8(6) 


0.5597 


27) 


[7,16] 


2.0(1.2) 


0.5483(28) 


0.02644(19) 










/? = 


2.10, Ksea — 


= 0.1374 










0.1357 


0.1357 


0.58900(50) 


[10,24] 


2.1(9) 


0.7281 


13) 


[11,24] 


1.6(8) 


0.8089(13) 


0.09906(47) 


0.1357 


0.1374 


0.51133(48) 


[10,24] 


1.6(8) 


0.6716 


17) 


[11,24] 


1.9(1.0) 


0.7614(19) 


0.07574(33) 


0.1367 


0.1367 


0.49686(47) 


[10,24] 


1.6(8) 


0.6612 


17) 


[11,24] 


1.8(9) 


0.7514(19) 


0.07174(31) 


0.1367 


0.1374 


0.46154(46) 


[10,24] 


1.4(8) 


0.6375 


20) 


[11,24] 


1.8(1.0) 


0.7240(22) 


0.06219(26) 


0.1374 


0.1374 


0.42401(46) 


[10,24] 


1.1(7) 


0.6133 


23) 


[11,24] 


1.7(1.0) 


0.6914(25) 


0.05267(22) 


0.1374 


0.1382 


0.37751(47) 


[10,24] 


0.9(6) 


0.5853 


27) 


[10,24] 


1.5(9) 


0.6450(29) 


0.04177(17) 


0.1374 


0.1385 


0.35879(48) 


[10,24] 


0.7(6) 


0.5757 


32) 


[10,24] 


1.5(9) 


0.6232(34) 


0.03766(16) 


0.1382 


0.1382 


0.32517(50) 


[10,24] 


0.7(6) 


0.5573 


36) 


[10,24] 


1.2(8) 


0.5835(37) 


0.03092(14) 


0.1385 


0.1385 


0.27972(59) 


[10,24] 


0.6(5) 


0.5367 


52) 


[10,24] 


1.1(7) 


0.5212(49) 


0.02271(13) 










/3 = 


2.10, /^sca — 


= 0.1382 










0.1357 


0.1357 


0.56682(64) 


[11,24] 


0.7(6) 


0.6923 


12) 


[10,24] 


0.6(5) 


0.8188(12) 


0.09396(53) 


0.1357 


0.1382 


0.44670(71) 


[10,24] 


1.2(7) 


0.6053 


16) 


[10,24] 


0.6(5) 


0.7379(18) 


0.05992(31) 


0.1367 


0.1367 


0.47282(67) 


[11,24] 


0.8(6) 


0.6227 


'13) 


[10,24] 


0.7(6) 


0.7593(15) 


0.06682(37) 


0.1367 


0.1382 


0.39191(75) 


[10,24] 


1.3(7) 


0.5689 


18) 


[10,24] 


0.9(7) 


0.6888(21) 


0.04640(26) 


0.1374 


0.1374 


0.39753(73) 


[10,24] 


1.3(7) 


0.5721 


16) 


[10,24] 


1.0(7) 


0.6949(20) 


0.04784(27) 


0.1374 


0.1382 


0.34943(78) 


[10,24] 


1.3(7) 


0.5417 


21) 


[9,24] 


1.1(8) 


0.6451(27) 


0.03693(23) 
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0.1382 


0.1382 


0.29459(85) 


[10,24] 


1.3 


7) 


0.5114(29) 


[9,24] 


1.1 


'J) 


0.5761(35) 


0.02613(18) 


0.1382 


0.1385 


0.27142(88) 


[10,24] 


1.2 


7) 


0.5003(34) 


[9,24] 


1.1 




0.5425(40) 


0.02208(17) 


0.1385 


0.1385 


0.24604(90) 


[10,24] 


1.2 


7) 


0.4887(42) 


[9,24] 


1.0 


;6) 


0.5034(46) 


0.01803(15) 










a 


= 2.20, Ksea = 0.1351 












0.1351 


0.1351 


0.49996(83) 


[11,24] 


0.5 


6) 


0.6260(23) 


[11,24] 


0.6 


;9) 


0.7987(28) 


0.08218(60) 


0.1351 


0.1358 


0.46683(86) 


[11,24] 


0.5 


6) 


0.6021(22) 


[11,24] 


0.4 


'J) 


0.7753(29) 


0.07232(52) 


0.1351 


0.1363 


0.44209(89) 


[11,24] 


0.5 


6) 


0.5849(22) 


[11,24] 


0.3 


[G) 


0.7558(30) 


0.06522(45) 


0.1351 


0.1368 


0.41631(96) 


[11,24] 


0.5 


6) 


0.5677(24) 


[11,24] 


0.4 


;5) 


0.7334(33) 


0.05806(38) 


0.1351 


0.1372 


0.3948(11) 


[11,24] 


0.5 


6) 


0.5538(30) 


[10,24] 


0.7 


;6) 


0.7130(41) 


0.05230(33) 


0.1358 


0.1358 


0.43202(89) 


[11,24] 


0.5 


6) 


0.5778(23) 


[11,24] 


0.4 


;6) 


0.7471(31) 


0.06247(43) 


0.1363 


0.1363 


0.37822(96) 


[11,24] 


0.5 


6) 


0.5428(29) 


[11,24] 


0.5 


;6) 


0.6969(39) 


0.04830(30) 


0.1368 


0.1368 


0.3174(11) 


[10.24] 


0.5 


7) 


0..-)087(;41) 


[10,21] 


1.3 


'1.0) 


0.6239(55) 


0.03413(28) 


0.1:i72 


0.i:i72 


0.2:)!)l)(i2) 


[i0.2i] 


0.9 


i.l) 


0. 1829(()()) 


[i0,2i] 




1.1) 


0.5:582(80) 


0.02271(27) 










/3 


= 2.20, Ksoa = 0.1358 












0.1351 


0.1351 


0.4879(15) 


[9,24] 


1.6 


1.1) 


0.6047(21) 


[10,24] 


0.6 


;7) 


0.8070(33) 


0.0797(10) 


0.1351 


0.1358 


0.4544(14) 


[9,24] 


1.4 


1.0) 


0.5803(26) 


[10,24] 


0.4 


;7) 


0.7830(38) 


0.06972(85) 


0.1358 


0.1358 


0.4190(13) 


[9,24] 


1.3 


1.0) 


0.5567(27) 


[9,24] 


0.6 


;i.5) 


0.7528(39) 


0.05971(77) 


0.1358 


0.1363 


0.3924(12) 


[9,24] 


1.3 


1.0) 


0.5391(30) 


[9,24] 


0.6 


;i.7) 


0.7278(41) 


0.05274(67) 


0.1358 


0.1368 


0.3643(12) 


[9,24] 


1.3 


1.0) 


0.5217(35) 


[9,24] 


0.7 


;i.4) 


0.6983(45) 


0.04561(61) 


0.1358 


0.1372 


0.3407(11) 


[9,24] 


1.3 


1.1) 


0.5084(42) 


[9,24] 


0.6 


;9) 


0.6701(53) 


0.03990(59) 


0.1363 


0.1363 


0.3643(12) 


[9,24] 


1.3 


1.0) 


0.5218(29) 


[9,19] 


0.4 


;5) 


0.6981(37) 


0.04544(55) 


0.1368 


0.1368 


0.3023(11) 


[9,24] 


1.4 


1.3) 


0.4855(33) 


[8,19] 


0.6 


;6) 


0.6227(40) 


0.03144(44) 


0.1372 


0.1372 


0.2434(11) 


[9,24] 


2.0 


1.5) 


0.4608(49) 


[8,19] 


0.6 


:6) 


0.5282(58) 


0.02004(40) 










/3 




.20, Ksea = 0.1363 












0.1351 


0.1351 


0.47893(93) 


[11,24] 


1.4 


1.11 


u.oy i^^zu ) 


[11,24] 


0.8 


'J) 


0.8099(22) 


0.07808(87) 


0.1351 


0.1363 


0.42030(98) 


[11,24] 


1.3 


1.1) 


0.5481(26) 


[11,24] 


1.4 


;i.o) 


0.7668(34) 


0.06115(71) 


0.1358 


0.1358 


0.41009(97) 


[11,24] 


1.3 


1.1) 


0.5406(26) 


[11,24] 


1.4 




0.7585(34) 


0.05839(67) 


0.1358 


0.1363 


0.3835(10) 


[11,24] 


1.3 


1.2) 


0.5223(31) 


[11,24] 


1.6 


;i.3) 


0.7341(42) 


0.05134(60) 


0.1363 


0.1363 


0.3554(10) 


[11,24] 


1.3 


1.2) 


0.5041(40) 


[10,24] 


1.6 


;i.6) 


0.7051(54) 


0.04431(54) 


0.1363 


0.1368 


0.3256(11) 


[11,24] 


1.3 


1.1) 


0.4862(49) 


[10,24] 


1.6 


;i.6) 


0.6698(66) 


0.03724(47) 


0.1363 


0.1372 


0.3002(12) 


[11,24] 


1.3 


1.1) 


0.4717(50) 


[10,24] 


1.4 


;i.3) 


0.6364(66) 


0.03156(39) 


0.1368 


0.1368 


0.2934(11) 


[11,24] 


1.4 


1.2) 


0.4676(54) 


[10,24] 


1.4 


;i.3) 


0.6274(71) 


0.03016(39) 


0.1372 


0.1372 


0.2338(13) 


[11,24] 


1.5 


1.4) 


0.4336(68) 


[10,24] 


1.6 




0.5392(81) 


0.01874(27) 










/3 




.20, Kaea = 0.1368 












0.1351 


0.1351 


0.4659(16) 


[10,24] 


4.3 


1.9) 


0.5715(24) 


[10,24] 


1.6 


;i.i) 


0.8152(45) 


0.0773(12) 


0.1351 


0.1368 


0.3810(24) 


[10,24] 


3.4 


1.8) 


0.5077(34) 


[9,23] 


0.8 


:8) 


0.7504(63) 


0.05265(88) 


0.1358 


0.1358 


0.3968(21) 


[10,24] 


4.0 


2.0) 


0.5181(31) 


[10,24] 


1.2 


;i.o) 


0.7658(61) 


0.05724(96) 


0.1358 


0.1368 


0.3419(23) 


[10,24] 


3.2 


1.8) 


0.4805(39) 


[9,23] 


0.7 


;8) 


0.7114(69) 


0.04278(69) 


0.1363 


0.1363 


3418(23^ 


10 24 


3.4 


1.9) 


0.4806(39) 


9 23 


0.7 


'8) 


7111(^70^ 


04289(701 


0.1363 


0.1368 


0.3115(23) 


[10,24] 


2.9 


1.8) 


0.4606(42) 


[9,23] 


0.7 


;8) 


0.6763(72) 


0.03569(56) 


0.1368 


0.1368 


0.2785(22) 


[10,24] 


2.5 


1.6) 


0.4407(44) 


[9,22] 


0.6 


:8) 


0.6320(70) 


0.02854(39) 


0.1368 


0.1372 


0.2496(23) 


[10,24] 


2.2 


1.4) 


0.4249(51) 


[9,22] 


0.8 


;i.o) 


0.5874(79) 


0.02284(39) 


0.1372 


0.1372 


0.2170(25) 


[10,24] 


2.0 


1.4) 


0.4086(66) 


[9,22] 


1.0 


;i.i) 


0.531(10) 


0.01721(31) 
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TABLE XXII. Baryon masses with degenerate quark combinations. 





mN 


[^miri) ^max] 








rriA 


[^min: ^max] 


XV^DF 








= 


1.80, Ksea = 


0.1409 








0.1409 


2.2596(49) 


[6,12] 




0.9(8) 




2.4151(71) 


[6,11] 


2.7(1.8) 


0.1430 


2.1085(62) 


[6,12] 




0.5(6) 




2.2853(94) 


[6,11] 


1.4(1.3) 


0.1445 


1.9950(82) 


[6,12] 




0.3(5) 




2.193(13) 


[6,11] 


0.6(8) 


0.1464 


1.8355(78) 


[5,11] 




0.3(5) 




2.061(13) 


[5,9] 


0.3(6) 


0.1474 


1.746(11) 


[5,11] 




0.4(6) 




1.999(19) 


[5,9] 


0.3(6) 








13 = 


1.80, Ksca = 


0.1430 








0.1409 


2.1797(53) 


[6,12] 




0.5(7) 




1.3326(88) 


[6,10] 


1.1(1.3) 


0.1430 


2.0237(70) 


[6,12] 




0.8(8) 




2.206(12) 


[6,10] 


0.9(1.1) 


0.1445 


1.9037(90) 


[6,12] 




1.3(1.1) 




2.104(16) 


[6,10] 


0.6(9) 


0.1464 


1.7397(93) 


[5,10] 




1.0(1.1) 




1.984(16) 


[5,9] 


0.2(5) 


0.1474 


1.644(11) 


[4,7] 




0.1(8) 




1.918(24) 


[5,9] 


0.4(8) 








/? = 


1.80, Ksoa = 


0.1445 








0.1409 


2.0963(56) 


[7,12] 




0.7(9) 




2.2364(58) 


[6,12] 


0.7(7) 


0.1430 


1.9362(70) 


[7,12] 




1.0(1.1) 




2.1033(78) 


[6,10] 


0.6(9) 


0.1445 


1.8126(93) 


[7,12] 




1.2(1.2) 




1.9978(67) 


[5,9] 


1.3(1.4) 


0.1464 


1.620(12) 


[6,9] 




1.7(1.9) 




1.871(12) 


[5,9] 


0.7(1.0) 


0.1474 


1.525(13) 


[5,8] 




1.9(1.7) 




1.826(24) 


[5,9] 


0.8(1.1) 








/3 = 


1.80, fvaea = 


0.1464 








0.1409 


1.9222(53) 


[6,12] 




1.6(1.1) 




2.0548(68) 


[6,12] 


0.7(8) 


0.1430 


1.7438(64) 


[6,12] 




1.4(1.0) 




1.9032(97) 


[6,12] 


1.0(9) 


0.1445 


1.6031(80) 


[6,12] 




0.9(8) 




1.789(14) 


[6,12] 


0.8(9) 


0.1464 


1.405(11) 


[5,9] 




0.1(3) 




1.655(12) 


[4,9] 


0.2(4) 


0.1474 


1.277(17) 


[4,9] 




0.3(5) 




1.572(21) 


[4,9] 


0.7(8) 








P = 


1.95, Ksea — 


0.1375 








0.1375 


1.7035(34) 


[8,14] 




0.4(5) 




1.8289(47) 


[7,14] 


0.6(7) 


0.1390 


1.5671(39) 


[8,14] 




0.4(6) 




1.7125(55) 


[7,14] 


0.5(6) 


0.1400 


1.4700(45) 


[8,14] 




0.5(7) 




1.6361(66) 


[7,14] 


0.5(6) 


0.1410 


1.3662(59) 


[8,14] 




0.5(7) 




1.5601(88) 


[7,13] 


0.6(7) 


0.1415 


1.3117(70) 


[8,14] 




0.4(6) 




1.523(11) 


[7,13] 


0.6(8) 








f3 = 


1.95, Ksca = 


0.1390 








0.1375 


1.6001(36) 


[8,13] 




0.4(7) 




1.7193(50) 


[8,16] 


0.4(5) 


0.1390 


1.4559(38) 


[7,13] 




0.5(6) 




1.6023(66) 


[8,16] 


0.5(5) 


0.1400 


1.3549(42) 


[7,13] 




0.6(7) 




1.5186(65) 


[7,14] 


0.3(5) 


0.1410 


1.2482(47) 


[7,12] 




0.7(8) 




1.4398(85) 


[7,14] 


0.2(4) 


0.1415 


1.1911(49) 


[6,12] 




0.6(7) 




1.401(11) 


[7,14] 


0.1(3) 








P = 


1.95, Ksea = 


0.1400 








0.1375 


1.5241(28) 


[6,15] 




1.0(7) 




1.6386(43) 


[6,15] 


1.2(9) 


0.1390 


1.3748(33) 


[6,15] 




1.4(9) 




1.5127(50) 


[6,15] 


1.1(8) 


0.1400 


1.2679(39) 


[6,15] 




2.1(1.2) 




1.4325(54) 


[6,15] 


1.0(8) 


0.1410 


1.1525(49) 


[6,15] 




2.4(1.3) 




1.3450(60) 


[6,11] 


0.6(8) 


0.1415 


1.0891(56) 


[6,15] 




2.2(1.1) 




1.3040(72) 


[6,11] 


0.7(9) 








13 = 


1.95, Ksea — 


0.1410 








0.1375 


1.4360(35) 


[6,16] 




1.9(1.0) 




1.5418(50) 


[6,16] 


1.8(1.0) 


0.1390 


1.2826(41) 


[6,15] 




1.8(1.0) 




1.4126(53) 


[6,14] 


1.5(1.0) 


0.1400 


1.1728(46) 


[6,15] 




1.4(9) 




1.3273(59) 


[6,14] 


1.1(8) 


0.1410 


1.0532(51) 


[6,13] 




1.2(1.0) 




1.2493(84) 


[6,12] 


0.7(8) 


0.1415 


0.9898(67) 


[6,13] 




1.7(1.1) 




1.206(11) 


[6,12] 


1.0(9) 








(3 = 


2.10, /^sca — 


0.1357 








0.1357 


1.1855(26) 


[11,20] 




2.2(1.4) 




1.2775(44) 


[11,22] 


0.8(7) 


0.1367 


1.0747(28) 


[11,20] 




1.8(1.2) 




1.1809(58) 


[11,22] 


0.7(6) 


0.1374 


0.9930(34) 


[11,20] 




1.5(1.0) 




1.1099(56) 


[10,19] 


1.1(8) 


0.1382 


0.8885(43) 


[10,20] 




1.5(1.0) 




1.0299(61) 


[8,18] 


0.9(7) 


0.1385 


0.8461(56) 


[10,20] 




1.4(9) 




1.0041(72) 


[8,13] 


1.2(1.1) 








13 = 


2.10, Ksea — 


0.1367 








0.1357 


1.1375(26) 


[11,23] 




1.9(1.0) 




1.2244(52) 


[12,22] 


2.5(1.2) 


0.1367 


1.0226(32) 


[11,23] 




1.3(9) 




1.1255(64) 


[11,22] 


2.6(1.2) 
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0.1374 


0.9363(35) 


[10,23] 




0.9(7) 




1.0562(71) 


[10,22] 


1.6(1.0) 


0.1382 


0.8311(45) 


[9,20] 




1.0(8) 




0.9770(85) 


[9,18] 


0.6(7) 


0.1385 


0.7888(55) 


[9,20] 




1.1(8) 




0.951(11) 


[9,16] 


0.4(6) 










2.10, /^sca — 


0.1374 








0.1357 


1.1046(34) 


[12,24] 




1.5(1.0) 




1.1797(44) 


[11,24] 


1.2(8) 


0.1367 


0.9873(33) 


[12,24] 




1.3(8) 




1.0781(60) 


[11,24] 


0.9(6) 


0.1374 


0.8955(35) 


[11,24] 




1.2(8) 




1.0089(59) 


[10,19] 


0.8(8) 


0.1382 


0.7866(49) 


[10,24] 




1.2(7) 




0.9301(93) 


[9,19] 


1.0(1.0) 


0.1385 


0.7438(67) 


[9,24] 




1.2(7) 




0.905(13) 


[9,19] 


1.3(1.1) 








f3 = 


2.10, /^sca — 


0.1382 








0.1357 


1.0526(33) 


[12,24] 




0.6(6) 




1.1233(41) 


[11,24] 


0.6(5) 


0.1367 


0.9319(36) 


[12,24] 




0.8(7) 




1.0168(55) 


[11,24] 


1.0(7) 


0.1374 


0.8383(38) 


[11,21] 




0.5(5) 




0.9389(59) 


[10,22] 


1.0(8) 


0.1382 


0.7204(42) 


[10,21] 




0.7(6) 




0.8887(92) 


[9,20] 


0.6(5) 


0.1385 


0.6680(65) 


[10,21] 




1.5(9) 




0.826(12) 


[9,19] 


0.7(6) 








P = 


2.20, /^sca — 


0.1351 








0.1351 


0.9330(76) 


[11,23] 




2.2(1.6) 




1.0219(65) 


[11,22] 


1.2(1.0) 


0.1358 


0.8463(84) 


[11,23] 




2.6(1.5) 




0.9474(73) 


[11,22] 


1.3(1.0) 


0.1363 


0.7843(79) 


[10,22] 




2.2(1.6) 




0.8931(74) 


[11,20] 


0.9(9) 


0.1368 


0.7155(90) 


[10,22] 




2.6(1.6) 




0.8444(89) 


[11,16] 


1.2(1.0) 


0.1372 


0.6540(88) 


[10,22] 




3.0(1.3) 




0.796(10) 


[10,15] 


1.3(9) 








13 = 


2.20, Ksca — 


0.1358 








0.1351 


0.9179(47) 


[12,23] 




1.2(1.1) 




0.9806(69) 


[11,24] 


1.2(1.1) 


0.1358 


0.8252(47) 


[11,23] 




1.5(1.1) 




0.9102(93) 


[11,23] 


1.1(1.0) 


0.1363 


0.7598(50) 


[10,22] 




0.8(8) 




0.8563(92) 


[11,21] 


0.6(8) 


0.1368 


0.6925(68) 


[10,22] 




1.1(1.1) 




0.8090(85) 


[9,21] 


0.7(8) 


0.1372 


0.6387(82) 


[10,20] 




1.2(1.0) 




0.782(12) 


[9,20] 


1.4(1.0) 








= 


2.20, Ksea — 


0.1363 








0.1351 


0.8915(41) 


[13,23] 




0.9(9) 




0.9610(46) 


[13,21] 


0.3(5) 


0.1358 


0.8028(43) 


[12,23] 




0.9(9) 




0.8882(46) 


[12,21] 


0.5(8) 


0.1363 


0.7349(42) 


[11,23] 




0.8(1.2) 




0.8343(62) 


[11,21] 


1.1(9) 


0.1368 


0.6663(60) 


[11,22] 




0.4(7) 




0.7775(97) 


[10,21] 


1.2(9) 


0.1372 


0.6014(98) 


[19,22] 




0.5(7) 




0.728(15) 


[10.21] 


0.9(8) 








J = 


2.20, K„.» = 


0.l:i()8 








0.1351 


860fi('36'l 


fl 1 241 




1.5(1.2) 




Q1 58Cfi4'l 

VJ . -L O Lj I Vjt: I 


[1 1 201 

[11, ^UJ 


2.1(1.2) 


0.1358 


0.7728(39) 


[10,24] 




1.2(1.0) 




0.8386(72) 


[10,20] 


1.6(1.0) 


0.1363 


0.7060(50) 


[10,24] 




1.1(8) 




0.7819(74) 


[9,20] 


1.4(9) 


0.1368 


0.6314(55) 


[9,23] 




0.7(6) 




0.7235(86) 


[9,19] 


1.5(1.0) 


0.1372 


0.5679(72) 


[9,23] 




0.6(9) 




0.676(13) 


[9,19] 


1.5(1.2) 
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TABLE XXIII. Baryon masses with non-degenerate quark combinations. 



'*val 


'*val 


'*val 




[tmin,tmax] X^/^DF 


rUA [tmin,tmax] X^/^DF 




[imin,tmax] / Nu¥ 










(3 = 1.80, 


Ksea = 0.1409 







0.1430 


0.1430 


0.1409 


2.1634(57) 


[6,12] 


0.6 


'J) 


2.1557(57) 


[6,12] 


0.6 


7) 


2.3280(84) 


[6,11] 


1.8 


1.4) 


0.1445 


0.1445 


0.1409 


2.0924(68) 


[6,12] 


0.4 


'6) 


2.0778(68) 


[6,12] 


0.5 


6) 


2.267(10) 


[6,11] 


1.1 


1.1) 


0.1464 


0.1464 


0.1409 


1.9981(64) 


rr 1 1 1 

[5,11] 


0.4 


6) 


1.9750(94) 


[6,11] 


0.3 


6) 


2.188(14) 


r/^ 1 1 1 

[6,11] 


0.4 


6) 


0.1474 


0.1474 


0.1409 


1 r\v r\-\ I cmW 

1.9501(80) 


[5,11] 


0.3 


5) 


1.919(14) 


[6,11] 


0.2 


r \ 

5) 


2.148(19) 


[6,11] 


0.3 


r \ 


0.1409 


0.1409 


0.1430 




[6,12] 


0.7 


'J) 


2.2135(53) 


[6,12] 


0.7 


8) 


2.3700(76) 


[^5 1 1 1 

[6,11] 


2.2 


1.6) 


0.1409 


0.1409 


0.1445 


2.1670(56) 


[6,12] 


0.6 


:7) 


2.1798(57) 


[6,12] 


0.6 


I) 


2.3392(83) 


[6,11] 


1.9 


1.5) 


0.1409 


0.1409 


0.1464 


2.1153(64) 


[6,12] 


0.5 


'6) 


2.1366(68) 


[6,12] 


0.7 


8) 


2.2993(94) 


[6,11] 


1.4 


1.2) 


0.1409 


0.1409 


0.1474 


2.0872(72) 


[6,12] 


0.4 


6) 


2.1144(77) 


[6,11] 


0.5 


7) 


2.277(11) 


r/^ 1 1 1 

[6,11] 


1.2 


1.1) 












P -- 


= 1.80, Ksoa = 0.1430 
















0.1409 


0.1409 


0.1430 


2.1246(58) 


[6,12] 


0.6 


8) 


2.1328(58) 


[6,12] 


0.6 


7) 


2.2893(97) 


[6,10] 


1.1 


1.2) 


0.1445 


0.1445 


0.1430 


■1 r\ A f\r\ ( ci A\ 

1.9490(84) 


[6,12] 


1.1 


'1.0) 


1.9402(82) 


[6,12] 


1.1 


1.0) 


r» io'?/i a\ 

2.137(14) 


[6,10] 


0.7 


1 A\ 

1.0) 


0.1464 


0.1464 


0.1430 


1.8505(77) 


[5,10] 


1.2 


1.2) 


1.8257(76) 


[5,10] 


0.8 


9) 


2.060(12) 


rr r\l 

[5,9] 


0.3 


6) 


0.1474 


0.1474 


0.1430 


1.7946(73) 


[4,10] 


1.0 


1.0) 


1.761(10) 


[5,10] 


1.1 


1.0) 


2.020(15) 


rr nl 

[5,9] 


0.3 


7) 


0.1430 


0.1430 


0.1409 


2.0813(64) 


[6,12] 


0.7 


'8) 


2.0721(64) 


[6,12] 


0.7 


8) 


2.246(11) 


[6,10] 


1.0 


1.2) 


0.1430 


0.1430 


0.1445 


1.9806(76) 


[6,12] 


0.9 


;9) 


1.9886(77) 


[6,12] 


1.0 


9) 


2.169(13) 


[6,10] 


0.8 


1.0) 


0.1430 


0.1430 


0.1464 


1.9230(64) 


[5,10] 


0.7 


'9) 


1.9426(91) 


[6,12] 


1.2 


1.0) 


2.1302(96) 


rr nl 

[5,9] 


0.5 


n\ 

9) 


0.1430 


0.1430 


0.1474 


1.8922(70) 


[5,10] 


0.9 


1.0) 


1.916(10) 


[6,12] 


1.1 


9) 


-1 -1 -I / -1 r\\ 

2.111(10) 


rc r»l 

[5,9] 


0.4 


8) 












P -- 


= 1.80, Ksea = 0.1445 
















0.1409 


0.1409 


0.1445 


1.9982(65) 


[7,12] 


0.6 


:8) 


2.0035(55) 


[6,12] 


2.3 


1.4) 


2.1595(68) 


[6,10] 


0.6 


9) 


0.1430 


0.1430 


0.1445 


1.8920(76) 


[7,12] 


0.9 


'1.1) 


1.8892(62) 


[6,12] 


2.2 


1 A\ 

1.4) 


2.0716(86) 


[6,10] 


0.7 


1.0) 


0.1464 


0.1464 


0.1445 


1.6892(65) 


rr 1^1 

[5,9] 


0.9 


1.2) 


1.6729(65) 


rr r»l 

[5,9] 


1.3 


1 a\ 

1.4) 


1.9103(99) 


rr ol 

[5,8] 


0.6 


1.2) 


0.1474 


0.1474 


0.1445 


1.6319(81) 


rr rkl 

[5,9] 


0.7 


1.0) 


1.6038(84) 


rr r»l 

[5,9] 


1.8 


1 a\ 

1.4) 


1.873(14) 


fr ol 

[5,8] 


0.6 


1.1) 


0.1445 


0.1445 


0.1409 


1.9103(62) 


[6,12] 


2.2 


1.3) 


1.8907(63) 


[6,12] 


2.0 


1 A\ 

1.4) 


2.0851(85) 


[6,10] 


0.7 


1 n\ 

1.0) 


0.1445 


0.1445 


0.1430 


1.8479(66) 


[6,12] 


2.1 


'1.3) 


1.8388(66) 


[6,12J 


2.1 


1 a\ 

1.4) 


2.0406(97) 


[6,10] 


0.8 


1.0) 


0.1445 


0.1445 


0.1464 


1.7321(59) 


rr 1^1 

[5,9] 


1.5 


1.4) 


1.7505(82) 


[6,12] 


1.7 


1 o\ 

1.3) 


1.9519(80) 


fr r\l 

[5,9] 


1.0 


1.2) 


0.1445 


0.1445 


0.1474 


1.6965(64) 


rr nl 

[5,9] 


1.1 


1.2) 


1.7242(94) 


[6,12] 


1.3 


1.1) 


1.9314(91) 


fr nl 

[5,9] 


1.0 


1.2) 












P -- 


= 1.80, Ksca = 0.1464 
















0.1409 


0.1409 


0.1464 


"1 n A c\ A f r*\ 

1.7494(66) 


[6,12] 


1.0 


;8) 


1.7765(65) 


r<" 111 

[6,11] 


1.8 


1 a\ 

1.4) 


1.9173(78) 


[5,10] 


0.2 


4) 


0.1430 


0.1430 


0.1464 


1.6253(78) 


[6,12] 


0.8 


'8) 


1.6465(79) 


[6,11] 


1.1 


1.0) 


1.8138(99) 


[5,10] 


0.3 


5) 


0.1445 


0.1445 


0.1464 


1.5320(97) 


[6,12] 


1.0 


9) 


1.5454(99) 


[6,11] 


0.8 


9) 


1.7379(91) 


r A nl 

[4,9] 


0.1 


3) 


0.1474 


0.1474 


0.1464 


-i 000/1 a\ 

1.333(14) 


r A n1 

[4,9] 


0.2 


.5) 


1.307(16) 


[5,10] 


0.5 


7) 


1.600(16) 


r A nl 

[4,9] 


0.5 


8) 


0.1464 


0.1464 


0.1409 


1.6086(77) 


rc 1 1 1 

[5,11] 


0.9 


;9) 


1.5700(92) 


[6,12] 


0.7 


7) 


■1 TOO / 1 n\ 

1.782(12) 


rc -1 nl 

[5,10] 


0.1 


4) 


0.1464 


0.1464 


0.1430 


1.5364(85) 


rr 1 1 1 

[5,11] 


0.8 


'9) 


1.509(11) 


[6,12] 


0.9 


9) 


1.7305(94) 


r A nl 

[4,9] 


0.1 


3) 


0.1464 


0.1464 


0.1445 


1.4809(92) 


rr 1 1 1 

[5,11] 


0.9 


1.0) 


1.465(13) 


[6,12] 


1.1 


1.1) 


1.694(10) 


r A nl 

[4,9] 


0.1 


4) 


0.1464 


0.1464 


0.1474 


1.3579(96) 


[4,9] 


0.2 


4) 


1.372(11) 


r,i r\l 

[4,9] 


0.1 


2) 


■1 f*C\ A / 1 0\ 

1.624(13) 


r^i nl 

[4,9] 


0.3 


6) 












P -- 


= 1.95, Ksoa = 0.1375 
















0.1390 


0.1390 


0.1375 


1.6161(39) 


[8,13] 


0.1 


'3) 


1.6103(36) 


[8,16] 


0.5 


6) 


1.7508(51) 


[7,14] 


0.5 


6) 


0.1400 


0.1400 


0.1375 


1.5571(41) 


[8,13] 


0.1 


'4) 


1.5437(39) 


[8,16] 


0.5 


5) 


1.7005(57) 


[7,14] 


0.5 


6) 


0.1410 


0.1410 


0.1375 


1.4971(47) 


[8,13] 


0.2 


'5) 


1.4727(46) 


[8,16] 


0.5 


6) 


1.6512(67) 


[7,14] 


0.5 


6) 


0.1415 


0.1415 


0.1375 


1.4669(53) 


[8,13] 


0.3 


5) 


1.4354(52) 


[8,16] 


0.6 


6) 


1.6271(74) 


[7,14] 


0.6 


7) 


0.1375 


0.1375 


0.1390 


1.6544(37) 


[8,13] 


0.1 


4) 


1.6611(37) 


[8,13] 


0.1 


3) 


1 TO 00 ( A C\\ 

1.7883(49) 


[7,14] 


0.5 


6) 


0.1375 


0.1375 


0.1400 


1.6212(38) 


[8,13] 


0.1 


'4) 


1.6335(39) 


[8,13] 


0.1 


3) 


1.7631(51) 


[7,14] 


0.5 


6) 


0.1375 


0.1375 


0.1410 


1.5868(39) 


[8,13] 


0.1 


4) 


1.6059(40) 


[8,13] 


0.2 


4) 


1 T'Ooo/roX 

1.7383(53) 


[7,14] 


0.4 


6) 


0.1375 


0.1375 


0.1415 


1 r ofw t \r\\ 


To 1 Ol 

[8.1.?] 


0.1 


4) 


1.. 5922(42) 


[8.13] 


0.2 


5) 


1 'i^^c / r r \ 
i.7zDz(55) 


[7,14] 


0.5 


6) 














= 1.9 


K = 0.i:i9() 
















0.1375 


0.1375 


0.1390 


1.5494(38) 


ft) 1 ol 

[8,13] 


0.5 


7) 


1.5566(38) 


[8,13] 


0.4 


7) 


1.6790(54) 


[8,16] 


0.4 


5) 


0.1400 


0.1400 


0.1390 


1.3928(41) 


rr- 1 ol 

[7,13] 


0.5 


6) 


1.3852(40) 


[7,13] 


0.6 


7) 


1.5452(62) 


rr^ 1 rl 

[7,15] 


0.4 


5) 


A 1 /1 1 n 
U.i4iU 


A 1 /1 1 n 
U.i4iU 


A 1 onn 
0.1390 


i.32oo(4o) 


[7,13] 


A C 

0.0 


6) 


1.3109(44) 


It 1 ol 


A ^ 


8) 


1 /1A0'7/''70\ 

i.4yj ( ( 1 1) 


\'V 1 c^l 

[7,15] 


A 

0.3 


0) 


0.1415 


0.1415 


0.1390 


1.2962(44) 


[6,13] 


0.5 


.6) 


1.2718(47) 


[7,12] 


0.8 


8) 


1.4688(81) 


[7,15] 


0.3 


4) 


0.1390 


0.1390 


0.1375 


1.5093(40) 


[8,13] 


0.5 


'J) 


1.5001(36) 


[7,13] 


0.5 


7) 


1.6391(59) 


[8,16] 


0.4 


5) 


0.1390 


0.1390 


0.1400 


1.4192(39) 


[7,13] 


0.6 


'7) 


1.4261(40) 


[7,13] 


0.5 


6) 


1.5715(59) 


[7,16] 


0.3 


4) 


0.1390 


0.1390 


0.1410 


1.3819(41) 


[7,12] 


0.6 


8) 


1.3966(42) 


[7,13] 


0.4 


6) 


1.5457(63) 


[7,16] 


0.3 


4) 


0.1390 


0.1390 


0.1415 


1.3623(42) 


[7,12] 


0.7 


9) 


1.3821(43) 


[7,13] 


0.4 


5) 


1.5335(66) 


[7,16] 


0.3 


4) 














= 1.95, Kse. = 0.1400 
















0.1375 


0.1375 


0.1400 


1.4350(31) 


[6,15] 


1.2 


;8) 


1.4480(30) 


[6,15] 


1.1 


8) 


1.5687(48) 


[6,15] 


1.2 


9) 


0.1390 


0.1390 


0.1400 


1.3363(35) 


[6,15] 


1.6 


;i.o) 


1.3431(34) 


[6,15] 


1.6 


1.0) 


1.4844(51) 


[6,15] 


1.1 


8) 



38 



0.1410 


0.1410 


0.1400 


1.1978(45) 


[6,15] 


2.5 


'1.3) 


1.1862 


'46) 


[6,15] 


2.3 


;i.3) 


1.3733(56) 


[6,11] 


0.5 


;8) 


0.1415 


0.1415 


0.1400 


1.1619(49) 


[6,15] 


2.5 


1.2) 


1.1416 


;51) 


[6,15] 


2.1 


;i.2) 


1.3467(61) 


[6,11] 


0.6 


;8) 


0.1400 


0.1400 


0.1375 


1.3646(34) 


[6,15] 


1.5 


1.0) 


1.3477 


;35) 


[6,15] 


1.5 


;9) 


1.4991(51) 


[6,15] 


1.1 


;8) 


0.1400 


0.1400 


0.1390 


1.3080(37) 


[6,15] 


1.9 




1.3002 


;37) 


[6,15] 


1.8 


;i.o) 


1.4567(52) 


[6,15] 


1.1 


;8) 


0.1400 


0.1400 


0.1410 


1.2254(42) 


[6,15] 


2.1 


'1.2) 


1.2355 


;4i) 


[6,15] 


2.4 


;i.3) 


1.4005(53) 


[6,11] 


0.6 


;8) 


0.1400 


0.1400 


0.1415 


1.2029(45) 


[6,15] 


2.0 


'1.2) 


1.2195 


;43) 


[6,15] 


2.5 


;i.3) 


1.3869(55) 


[6,11] 


0.6 


;8) 












P-- 


— 1.95, /^sea — 


= 0.1410 














0.1375 


0.1375 


0.1410 


1.3059(42) 


[6,15] 


1.8 


'1.0) 


1.3260 


;4i) 


[6,14] 


1.4 


;i.o) 


1.4441(52) 


[6,14] 


1.4 


;9) 


0.1390 


0.1390 


0.1410 


1.2018(45) 


[6,15] 


1.5 


;9) 


1.2167 


'45) 


[6,14] 


1.0 


;8) 


1.3574(57) 


[6,14] 


1.2 


;9) 


0.1400 


0.1400 


0.1410 


1.1290(48) 


[6,15] 


1.4 


'9) 


1.1381 


;46) 


[6,13] 


0.7 


;8) 


1.3003(65) 


[6,14] 


0.9 


;8) 


0.1415 


0.1415 


0.1410 


1.0156(60) 


[6,13] 


1.5 


1.0) 


1.0072 


;62) 


[6,13] 


1.6 




1.2191(98) 


[6,12] 


1.0 


;9) 


0.1410 


0.1410 


0.1375 


1.2050(48) 


[6,14] 


1.0 


.8) 


1.1722 


;44) 


[6,13] 


0.7 


'J) 


1.3464(62) 


[6,14] 


0.8 


;8) 


0.1410 


0.1410 


0.1390 


1.1440(49) 


[6,14] 


1.0 


:8) 


1.1217 


;45) 


[6,13] 


0.8 




1.3026(66) 


[6,14] 


0.8 


[7) 


0.1410 


0.1410 


0.1400 


1.0997(48) 


[6,13] 


0.8 


'8) 


1.0874 


'47) 


[6,13] 


1.0 


;9) 


1.2737(72) 


[6,14] 


0.8 


;8) 


0.1410 


0.1410 


0.1415 


1.0286(56) 


[6,13] 


1.5 


1.1) 


1.0364 


;55) 


[6,13] 


1.3 


;i.o) 


1.2327(91) 


[6,14] 


1.2 


;9) 



P = 2.10, Ksoa = 0.1357 



0.1367 
0.1374 
0.1382 
0.1385 
0.1357 
0.1357 
0.1357 
0.1357 



0.1367 
0.1374 
0.1382 
0.1385 
0.1357 
0.1357 
0.1357 
0.1357 



0.1357 
0.1357 
0.1357 
0.1357 
0.1367 
0.1374 
0.1382 
0.1385 



1.1157(28) 
1.0661(32) 
1.0059(34) 
0.9845(41) 
1.1460(26) 
1.1173(26) 
1.0832(28) 
1.0701(28) 



11,20] 
11,20] 
10,20] 
10,20] 
11,20] 
11,20] 
11,20] 
11,20] 



2.2 
2.2 
2.2 
2.3 
1.9 
1.7 
1.6' 
1.6 



1.3) 
1.2) 
1.2) 
1.2) 
1.3) 
1.2) 
1.1) 
1.0) 



1.1086 
1.0513 
0.9808 
0.9523 
1.1521 
1.1285 
1.1013 
1.0912 



27) 
28) 
33) 
37) 
27) 
28) 
33) 
36) 



11,20] 
11,20] 
11,20] 
11,20] 
11,20] 
11,20] 
11,20] 
11,20] 



1.8 
1.5 
1.1 
1.0 
2.3 
2.4 
2.3 
2.2 



1.2) 
1.0) 
0.9) 
0.9) 
1.4) 
1.4) 
1.4) 
1.3) 



1.2098(41) 
1.1658(47) 
1.1125(50) 
1.0952(58) 
1.2409(38) 
1.2190(40) 
1.1940(45) 
1.1853(49) 



[10,19] 
[10,19] 
[9,18] 
[9,18] 
[10,19] 
[10,19] 
[10,18] 
[10,18] 



l.Oi 
l.li 
0.9i 
0.8i 
l.li 
l.Oi 
0.8i 
0.7i 



2.10, Ks. 



0.1367 



0.1357 

0.1374 
0.1382 
0.1385 
0.1367 
0.1367 
0.1367 
0.1367 



0.1357 

0.1374 
0.1382 
0.1385 
0.1367 
0.1367 
0.1367 
0.1367 



0.1367 

0.1367 
0.1367 
0.1367 
0.1357 
0.1374 
0.1382 
0.1385 



1.0970(27) 

0.9683(33) 
0.9061(40) 
0.8832(44) 
1.0647(30) 
0.9909(31) 
0.9542(34) 
0.9402(36) 



11,23] 

10,23] 
[9,20] 
[9,20] 
11,23] 
10,23] 
10,23] 
10,23] 



1.6 

1.0 
1.2 
1.3 
1.6 
1.1 
0.8 
0.7 



1.0) 

7) 
8) 
8) 
9) 
8) 
7) 
6) 



1.1026 

0.9627 
0.8890 
0.8599 
1.0584 
0.9972 
0.9680 
0.9578 



28) 

37) 
39) 

41) 
29) 
36) 

34) 
37) 



11,23] 

11,23] 
10,21] 
10,19] 
11,23] 
11,23] 
10,23] 
10,23] 



1.8 

1.0 
0.7 
0.8 
1.4 
1.3 
1.2 
1.2 



1.0) 

8) 
7 
7 
9 



1.1929(57) 

1.0788(66) 
1.0266(70) 
1.0087(80) 
1.1621(63) 
1.1007(61) 
1.0770(69) 
1.0691(74) 



[12,22] 

[10,22] 
[9,21] 
[9,21] 
[12,22] 
[10,22] 
[10,22] 
[10,22] 



2.3i 
1.9i 
1.2i 
l.Oi 
2.1i 
2.1i 
1.8i 
1.7i 



P 



2.10, Ks. 



0.1374 



0.1357 


0.1357 


0.1374 


1.0339(36) 


[12,24] 


1.5 


'1.0) 


1.0418 


;33) 


12,23] 


0.8 


8) 


1.1224(54) 


[11,24] 


0.9 


;6) 


0.1367 


0.1367 


0.1374 


0.9530(30) 


[11,24] 


1.8 


'1.0) 


0.9576 


;30) 


11,23] 


1.2 


8) 


1.0531(50) 


[10,19] 


0.8 


;8) 


0.1382 


0.1382 


0.1374 


0.8303(44) 


[10,24] 


1.0 


7) 


0.8177 


;38) 


10,23] 


1.3 


8) 


0.9549(76) 


[9,19] 


0.9 


;i.o) 


0.1385 


0.1385 


0.1374 


0.8065(48) 


[9,20] 


1.1 


8) 


0.7839 


;49) 


10,23] 


1.3 


7) 


0.9395(94) 


[9,19] 


1.0 


;i-i) 


0.1374 


0.1374 


0.1357 


0.9728(33) 


[11,24] 


1.0 


'7) 


0.9601 


;32) 


11,23] 


2.0 


1.2) 


1.0641(50) 


[10,19] 


0.9 


;8) 


0.1374 


0.1374 


0.1367 


0.9286(33) 


[11,24] 


1.1 


'7) 


0.9226 


;3i) 


11,23] 


1.7 


1.0) 


1.0284(53) 


[9,19] 


0.9 


;7) 


0.1374 


0.1374 


0.1382 


0.8553(32) 


[10,24] 


1.3 


8) 


0.8648 


;36) 


10,23] 


0.9 


6) 


0.9789(64) 


[9,19] 


0.8 


;8) 


0.1374 


0.1374 


0.1385 


0.8395(35) 


[10,24] 


1.4 


8) 


0.8540 


;42) 


10,23] 


1.0 


7) 


0.9705(70) 


[9,19] 


0.9 


;9) 












P-- 


— 2.10, Ksea ~ 


= 0.1382 














0.1357 


0.1357 


0.1382 


0.9425(35) 


[11,24] 


0.8 


'7) 


0.9568 


;33) 


10,24] 


0.9 


6) 


1.0352(59) 


[11,24] 


0.9 


:6) 


0.1367 


0.1367 


0.1382 


0.8596(38) 


[11,24] 


0.9 


'7) 


0.8717 


'35) 


10,24] 


0.9 


6) 


0.9614(56) 


[10,22] 


1.0 


;7) 


0.1374 


0.1374 


0.1382 


0.7969(33) 


[10,21] 


0.6 


6) 


0.8058 


;35) 


10,21] 


0.4 


5) 


0.9098(70) 


[10,22] 


0.9 


;7) 


0.1385 


0.1385 


0.1382 


0.6894(55) 


[10,21] 


1.1 


8) 


0.6824 


;48) 


[9,21] 


1.0 


7) 


0.836(11) 


[9,20] 


0.7 


;6) 


0.1382 


0.1382 


0.1357 


0.8529(37) 


[10,21] 


0.3 


5) 


0.8270 


;37) 


10,24] 


0.7 


6) 


0.9440(64) 


[10,22] 


0.9 


;7) 


0.1382 


0.1382 


0.1367 


0.8047(37) 


[10,21] 


0.4 


5) 


0.7863 


;37) 


10,24] 


0.7 


7) 


0.9073(75) 


[10,22] 


0.8 


'J) 


0.1382 


0.1382 


0.1374 


0.7677(38) 


[10,21] 


0.5 


5) 


0.7561 


(35) 


10,21] 


0.6 


6) 


0.8814(74) 


[9,22] 


0.8 


'J) 


0.1382 


0.1382 


0.1385 


0.7005(42) 


[9,21] 


0.8 


.6) 


0.7066 


;42) 


[9,21] 


0.7 


6) 


0.8448(99) 


[9,20] 


0.7 


;6) 












P-- 


— 2.20, Ksea ~ 


= 0.1351 














0.1358 


0.1358 


0.1351 


0.8784(78) 


[11,23] 


2.3 


'1.4) 


0.8729 


:81) 


10,23] 


2.3 


1.4) 


0.9723(70) 


[11,22] 


1.3 


;9) 


0.1363 


0.1363 


0.1351 


0.8382(85) 


[11,23] 


2.5 


1.5) 


0.8279 


;93) 


10,23] 


3.0 


1.6) 


0.9372(78) 


[11,22] 


1.5 


;i.o) 


0.1368 


0.1368 


0.1351 


0.801(11) 


[10,23] 


2.6 


1.7) 


0.7834 


;8i) 


10,23] 


3.5 


1.8) 


0.9016(81) 


[11,19] 


1.1 


;i.o) 


0.1372 


0.1372 


0.1351 


0.768(11) 


[10,23] 


2.3 


,1.6) 


0.7452 


;72) 


10,23] 


2.6 


1.4) 


0.874(10) 


[11,19] 


1.5 


;i.o) 


0.1351 


0.1351 


0.1358 


0.9022(82) 


[11,23] 


2.3 


;i.5) 


0.9064 


;72) 


10,23] 


2.0 


1.3) 


0.9963(67) 


[11,22] 


1.2 


:9) 


0.1351 


0.1351 


0.1363 


0.8793(86) 


[11,23] 


2.6 


'1.6) 


0.8880 


'73) 


10,23] 


2.0 


1.3) 


0.9792(70) 


[11,22] 


1.3 


;9) 


0.1351 


0.1351 


0.1368 


0.8562(91) 


[10,23] 


3.0 


.1-8) 


0.8698 


;78) 


10,23] 


1.9 


1.4) 


0.9623(74) 


[11,21] 


1.3 


;8) 
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0.1351 0.1351 0.1372 0.8397(91) [10,23] 3.8.(2.2) 0.8561(85) [10,23] 1.9(1.4) 0.9489(80) [11,21] 1.5(1.0) 

/3 = 2.20, Ksoa = 0.1358 



0.1351 


0.1351 


0.1358 


0.8874(48) 


[12,23] 


1.3 


1.1) 


0.8859(44) 


[11,23] 


1.6 




0.9565(77) 


[11,24] 


1.2 




0.1363 


0.1363 


0.1358 


0.7842(46) 


[10,22] 


0.8 


8) 


0.7795(46) 


[10,22] 


1.0 


;8) 


0.8748(84) 


[10,21] 


0.5 


;7) 


0.1368 


0.1368 


0.1358 


0.7438(56) 


[10,22] 


0.7 


8) 


0.7321(55) 


[10,22] 


1.1 


;i.o) 


0.8443(84) 


[10,21] 


0.6 


;8) 


0.1372 


0.1372 


0.1358 


0.7144(70) 


[10,22] 


0.9 


9) 


0.6930(57) 


[10,22] 


1.0 


;i.2) 


0.8224(82) 


[9,20] 


0.9 


;9) 


0.1358 


0.1358 


0.1351 


0.8616(49) 


[12,22] 


1.0 


9) 


0.8523(44) 


[11,23] 


1.7 


;i.2) 


0.9326(87) 


[11,24] 


1.1 




0.1358 


0.1358 


0.1363 


0.8016(43) 


[10,22] 


1.1 


8) 


0.8057(44) 


[10,22] 


0.9 


;8) 


0.8914(83) 


[10,21] 


0.6 


'J) 


0.1358 


0.1358 


0.1368 


0.7776(45) 


[10,22] 


1.0 


9) 


0.7863(48) 


[10,22] 


0.7 


;7) 


0.8764(84) 


[10,21] 


0.6 


;8) 


0.1358 


0.1358 


0.1372 


0.7586(48) 


[10,22] 


1.1 


1.0) 


0.7719(57) 


[10,22] 


0.6 


;8) 


0.8664(80) 


[9,21] 


0.8 


;8) 












!3-- 


= 2.20, Ksoa = 0.1363 














0.1351 


0.1351 


0.1363 


0.8366(44) 


[12,23] 


0.9 


9) 


0.8446(43) 


[12,23] 


0.8 


;8) 


0.9203(50) 


[12,20] 


0.4 


;7) 


0.1358 


0.1358 


0.1363 


0.7792(45) 


[12,23] 


0.8 


9) 


0.7831(45) 


[12,23] 


0.8 


;9) 


0.8702(56) 


[12,20] 


0.8 


;9) 


0.1368 


0.1368 


0.1363 


0.6939(54) 


[11,22] 


0.4 


6) 


0.6865(52) 


[10,22] 


0.6 


;8) 


0.7960(97) 


[11,19] 


1.2 


;i.2) 


0.1372 


0.1372 


0.1363 


0.6595(70) 


[11,22] 


0.2 


4) 


0.6400(63) 


[10,22] 


0.5 


;8) 


0.764(12) 


[10,19] 


1.1 




0.1363 


0.1363 


0.1351 


0.7940(46) 


[12,23] 


0.8 


9) 


0.7847(46) 


[12,23] 


0.8 


;9) 


0.8778(56) 


[12,20] 


0.8 


;9) 


0.1363 


0.1363 


0.1358 


0.7612(46) 


[12,23] 


0.8 


1.0) 


0.7570(48) 


[12,23] 


0.8 


;i.o) 


0.8528(60) 


[12,20] 


1.0 


;i.o) 


0.1363 


0.1363 


0.1368 


0.7107(45) 


[11,22] 


0.6 


8) 


0.7166(47) 


[10,22] 


0.6 


;7) 


0.8143(73) 


[11,20] 


1.3 


;i.i) 


0.1363 


0.1363 


0.1372 


0.6883(52) 


[11,22] 


0.5 


8) 


0.7003(55) 


[10,22] 


0.5 


;6) 


0.8002(85) 


[10,20] 


1.3 


;i.o) 












/3 = 


= 2.20, Ksoa = 0.1368 














0.1351 


0.1351 


0.1368 


0.7827(38) 


[10,24] 


1.2 


9) 


0.7956(41) 


[10,24] 


1.6 


;i.2) 


0.8536(69) 


[10,20] 


1.4 


;9) 


0.1358 


0.1358 


0.1368 


0.7242(46) 


[10,24] 


1.0 


7) 


0.7329(48) 


[10,24] 


1.3 


;i.o) 


0.8011(74) 


[9,20] 


1.3 


;8) 


0.1363 


0.1363 


0.1368 


0.6793(51) 


[10,23] 


0.9 


7) 


0.6847(52) 


[10,24] 


1.1 


;8) 


0.7624(74) 


[9,20] 


1.4 


;9) 


0.1372 


0.1372 


0.1368 


0.5935(64) 


[9,22] 


0.4 


5) 


0.5858(64) 


[9,23] 


0.5 


;7) 


0.691(12) 


[9,19] 


1.5 




0.1368 


0.1368 


0.1351 


0.7204(53) 


[10,23] 


1.2 


8) 


0.7060(46) 


[10,24] 


0.8 


;7) 


0.7898(74) 


[9,20] 


1.1 


;8) 


0.1368 


0.1368 


0.1358 


0.6858(53) 


[10,23] 


1.0 


7) 


0.6767(51) 


[10,24] 


0.7 


:6) 


0.7625(75) 


[9,20] 


1.2 


;9) 


0.1368 


0.1368 


0.1363 


0.6599(54) 


[10,23] 


0.9 


7) 


0.6541(47) 


[9,24] 


0.7 


;6) 


0.7424(78) 


[9,20] 


1.3 


;i.o) 


0.1368 


0.1368 


0.1372 


0.6076(57) 


[9,23] 


0.5 


6) 


0.6138(64) 


[9,23] 


0.7 


;7) 


0.707(10) 


[9,20] 


1.3 


;i.o) 



TABLE XXIV. 1x1 and 1x2 Wilson loops in full QCD at each simulated sea quark mass and extrapolated to the chiral 
limit. 



K 



1.80 



0.1409 
0.1430 
0.1445 
0.1464 
mps = 



0.490527(30) 
0.495049(39) 
0.499361(37) 
0.507204(57) 
0.51471(34) 



0.232159(35) 
0.237880(53) 
0.243370(49) 
0.253308(78) 
0.26274(45) 



1.95 



0.1375 
0.1390 
0.1400 
0.1410 
mps = 



0.553355(20) 
0.556667(21) 
0.559143(21) 
0.561884(27) 
0.56518(20) 



0.305089(27) 
0.309890(31) 
0.313473(34) 
0.317457(36) 
0.32228(31) 



2.10 



0.1357 
0.1367 
0.1374 
0.1382 

mps = 



0.5980283(76) 
0.5992023(76) 
0.6000552(67) 
0.6010819(84) 
0.602197(64) 



0.362139(12) 
0.363979(12) 
0.365297(10) 
0.366883(13) 
0.36862(10) 



2.20 



0.1351 
0.1358 
0.1363 

0.1368 

mps = 



0.620027(10) 
0.620616(7) 
0.621035(8) 

0.621490(8) 
0.62233(22) 



0.390976(16) 
0.391911(12) 
0.392570(11) 

0.393289(12) 
0.39465(26) 
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TABLE XXV. 1 X 1 and 1 X 2 Wilson loops in quenched QCD. 



p 






2.187 


0.5921968(62) 


0.3438930(92) 


2.214 


0.5991994(56) 


0.3533512(86) 


2.247 


0.6072343(59) 


0.3642759(92) 


2.281 


0.6149775(51) 


0.3748875(80) 


2.334 


0.6261248(48) 


0.3902635(80) 


2.416 


0.6415604(20) 


0.4117512(34) 


2.456 


0.6484512(17) 


0.4214117(29) 


2.487 


0.6535537(21) 


0.4286002(36) 


2.528 


0.6600072(20) 


0.4377232(34) 


2.575 


0.6670422(18) 


0.4477145(25) 



TABLE XXVL Numerical values for coupling constants, Z factors and improvement coefficients in full QCD. Z factors and 
improvement coefficionts are evaluated using gL—il/a). We also quote run factors used for running quark masses from fi = 1/a 
to 11 = 2 GeV with the three-loop beta function. 













Zm 




Zp 


hp 


Za 


Za/Zp 


hA 


CA 


Zv 


hv 


run factor 


1.80 


3.155 


2.185 
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FIG. 1. Schematic plot for the choice of sea and valence hopping parameters. For circles at the points Kvai = '^soa the 
corresponding pseudoscalar to vector meson mass ratio is indicated. 
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FIG. 2. Effective mass plots for psoudoscalar, vector, nuclcon and A clianncls witli degenerate valence hopping parameters 
Hvai = Usea. = 0.1445 at /3 = 1.8. Circles represent effective masses obtained when all quark propagators are calculated with 
point sources. For squares all quark propagators have smeared sources and triangles are for mixed combinations of sources. 
Solid lines denote the results from correlated mass fits to smeared source hadron propagators. Dashed lines show the one 
standard deviation error band determined by jack-knife analysis with a bin size of 10 configurations. 
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FIG. 3. Effective mass plots for pseudoscalar, vector, nucleon and A channels with degenerate valence hopping parameters 
Kvai = Kaea = 0.1400 at /? = 1.95. Symbols have the same meaning as in Fig. 0. 
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FIG. 4. Effective mass plots for pseudoscalar, vector, nucleon and A channels with degenerate valence hopping parameters 
= 0.1374 at /? = 2.1. Symbols have the same meaning as in Fig. |[ 
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FIG. 5. Effective mass plots of potential data at r = L/4 for sea quark mass corresponding to mps/my ~ 0.7. The scale is 
fixed from p meson mass at the physical point. 




FIG. 6. Static quark potential on 24 x 48 lattice at Ksea = 0.1374. Both vertical and horizontal lines are normalized by the 
Sommer scale ro. The solid line represents the fit curve of Eq. (Ba). 
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FIG. 7. Two examples of autocorrelation function (lower symbols) and cumulative autocorrelation time (upper symbols) for 
Afinv Errors are determined with the jack-knife method. Also plotted are two lines y{t) = t/4 and y{t) = t/10 within which a 
plateau of r'^"™(/;) can be observed. 
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FIG. 8. Relative errors of the pseudoscalar meson mass as a function of the bin size. Two examples are shown, each at the 
lightest sea quark mass of mps/mv »i 0.6. Data at the heaviest valence quark mass are represented by filled symbols and the 
ones from the lightest valence quark mass with open symbols. 
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FIG. 13. Chiral extrapolations ol octet baryon masses. Plots only show E-like octet baryons. Lines are from fits with 
Eq. dHi). 
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FIG. 14. Chiral extrapolations of string tension and Sommer scale. 
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FIG. 15. Partially quenched spectrum at the physical sea quark mass. Lines are obtained from Eqs. ( p8[ ) 
fixing /isoa = rn^. The strange spectrum, marked with symbols on the lines, is obtained using A^k as input. 
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FIG. 19. Chiral extrapolations ol vector meson masses in quenched QCD. Lines are Irom linear fits with Eq. (p 



1.7 



1.5 



1.3 



> 1.1 



0.9 



0.7 



0.5 



0^ ' ' 



0.0 



0.2 



0.4 

mp23(Kv,Kv) 



0.6 



FIG. 20. Chiral extrapolations of octet baryon masses in quenched QCD. While fits have been made to E and A type baryons 
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FIG. 22. Meson masses in quenched QCD with improved and standard actions. 
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FIG. 23. Octet baryon masses in quenched QCD with improved and standard actions. The strange spectrum is determined 
with K input (left figure) or input (right figure). 



56 



1.8 



1.6 



1.2 



no 



K input 



5§f 



> 
o 



^14- 



A O 



0.0 



0.5 

a [GeV-'] 



□ Standard Action 
o Improved Action 

1.0 



1.8 



1.6 



> 

O 



1.2 



no 



1 4 - 



A O 



0.0 



4) input 







.55. 



0.5 

a [GeV-'] 



□ Standard Action 
o Improved Action 

1.0 



FIG. 24. Dccuplct baryon masses in quenched QCD with improved and standard actions. The strange spectrum is determined 
with K input (left figure) or (j) input (right figure). 
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FIG. 29. Average up and down quark mass for three different definitions in full QCD. Lines are from linear extrapolations 
to the continuum limit made separately for each definition. 
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FIG. 30. Strange quark mass for three different definitions and two different experimental inputs in full QCD. Lines are from 
linear extrapolations to the continuum limit made separately for each definition. 
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FIG. 31. Comparison of average up and down quark mass in quenched and full QCD. Lines are from combined linear 
continuum extrapolations as described in the text. 
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FIG. 35. Lattice spacing dependence of pseudoscaJar decay constants /,r and /k in full QCD (filled circles) and quenched 
QCD with improved actions (large open circles) or standard action (small open squares). The strange quark mass used in the 
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FIG. 36. Comparison of /k/Ztt — 1 in full and quenched QCD. Fit lines are linear for all data. 
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FIG. 39. Lattice spacing dependence of decay constants Fp and in full and quenched QCD. 
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FIG. 40. Ratio of pseudoscalar and vector decay constants in full and quenched QCD. 
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FIG. 41. Perturbative and non-perturbative Z-factors for vector current at zero quark mass. 
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